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Fracture  failure  In  multi-layer  epoxy-baaed  composite  laminates  seldom 
begins  with  breaking  of  the  load-carrying  reinforcing  fibers.  Rather,  nail 
cracks  In  the  matrix-phase  of  the  material  appear  flrat,  usually  at  a  low 
stress  level.  Transverse  cracks  and  free-edge  da laminations  are  among  the 
most  often  observed  siib-laainate  fracture  eodes  in,  e.g. ,  glass-epoxy  and  graph¬ 
ite-epoxy  composite  laminates.  Although  these  matrix-dominated  crackings  do 
not  generally  cause  catastrophic  laminate  failure,  these  are  capable  of  propa¬ 
gating  Into  large  proportions  under  sustained  loading.  Thus,  the  growth  be¬ 
havior  of  the  sub-laminate  cracks  has  been  the  subject  of  Intensive  research 
In  recent  years. 

From  a  fracture  analysis  point  of  view,  the  growth  mechanisms  of  sub- 
laminate  crackings  may  be  described  at  two  different  dimensional  levels.  The 
first  approach  Is  at  the  micromachanlcs  level.  At  this  level,  quantities  of 
dimension  comparable  to  the  fiber  diameter  must  be  considered.  Thus,  the 
analysis  would  distinguish  the  exact  physical  configurations  in  the  fiber-matrix 
micro structure;  and  as  far  as  fracture  mechanisms  are  concerned,  the  analysis 
would  also  include  the  effects  of  micro-defects,  such  as  voids,  fiber-matrix 
disbonds,  etc.  Since  the  latter  may  distribute  randomly  throughout  the 
material,  any  fractura  analysis  method  based  on  micromechanics  must  also  address 
the  statistical  characteristics  of  the  material  micro-defects. 

Owing  to  the  extrema  complexity  In  die  micros tructures  of  any  given  fibrous 
composite,  fracture  analysis  of  sub- laminate  cracks  based  on  micromechanies 
Is  almost  Impossible  if  not  impractical.  Nevertheless,  the  approach  has  been 
applied  extensively  to  study  failure  mechanisms  of  unidirectional  fibrous 
composites,  where  fiber  failure  is  predominant. 


At  a  much  larger  dimensional  level,  a  phenomenological  fracture  analysis 
may  be  performed.  In  this  approach,  each  of  the  material  layers  In  the  lami¬ 
nate  Is  Idealised  as  a  homogeneous  anisotropic  medium,  where  the  fiber-matrix 
distinction  disappears .  The  physical  properties  of  the  Individually  Idealised 
layers  are  represented  by  some  average  values  from  that  of  the  fiber  and  the 
matrix  phases.  Consequently,  all  of  the  micro-defects  in  the  idealised  layer 
become  indistinguishable;  and  their  resultant  macroscopic  effect  on  failure  may 
be  regarded  aa  some  Inherent  material  property: 

Generally,  a  crack  of  dimension  much  larger  than  the  fiber  diameter  will 
be  recognized  at  the  macroscopic  level,  such  as  a  transverse  crack  In  a  layer, 
or  a  delamination  crack  In  a  layer-interface.  At  this  level  of  analysis, 
methods  of  the  classical  fracture  mechanics  may  be  used,  where  the  effective 
properties  of  the  material  layers  are  required  in  order  to  perform  the  necessary 
stress  analysis.  In  fact.  It  Is  at  this  level,  material  characterisations  of 
various  composite  material  systems  are  conducted  In  practice. 

Research  on  the  sub-laminate  fracture  growth  in  laminates  has  been 
mostly  experimental.  And  serious  analytical  studies  of  the  various  growth 
mechanisms  from  either  a  microscopic  approach  or  a  macroscopic  approach  have 
been  very  limited  (see,  e.g.,  the  review  on  fracture  analysis  of  composites 
given  recently  In  a  treatise  by  Cherepanov  [1]). 

In  what  follows,  a  more  detailed  discussion  will  be  given  to  the  two 
types  of  sub- laminate  crackings,  which  have  been  the  focus  of  considerable 
recent  Interest.  Specifically  we  shall  discuss  the  growth  behavior  of 
transverse  cracks  and  free  edge  delamination  in  graphite -epoxy  composite 
laminates. 
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ftwmm  cracking  is  pa  wily  mtdttd  with  tht  ln-plsae  twill 
•trap  which  la  normal  to  tiw  fibers  Is  •  gtp  material  layer.  Consular,  for 
instance,  the  l0/90]#  laminate  under  uniaxial  tension  as  shown  in  Tig.  1.1. 
Multiple  transverse  cracks  are  shown  to  occur  in  the  90*-layer.  is suns  that 
the  90*-layer  has  a  uniform  tensile  strength,  then  multiple  cracks  would 
appear  simultaneously  at  a  regular  spacing  when  some  critical  stress  is  reached 
transverse  to  the  fibers.  The  crack  spacing  is  determined  by  the  shear  stress 
transfer  zone  on  the  90/0  interface  where  the  transverse  crack  terminates. 

Thus,  first  intuition  would  suggest  that  the  critical  stress  in  the  90*-layer 
at  the  onset  of  multiple  cracking  must  be  the  tensile  strength  of  the  90*- 
layer.  The  latter  is,  of  course,  a  material  property  of  the  90* -layer. _ 


P 


•igure  1.*  Schematics  of  Multiple  Transverse  Cracks 
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Actually,  experiments  have  Indicated  that  the  tensile  stress  in  the 
90*-layer  at  onset  of  transverse  cracking  varies  greatly  with  the  thickness 
of  the  90 ‘-layer  itself.  The  general  trend  is  that  the  critical  tensile  stress 
in  90*-layer  increases  as  the  actual  thickness  of  the  90* -layer  decreases. 

This  layer-thickness  effect  was  first  documented  experimentally  by  Bader, 
et  al.  [2],  who  tested  lsmlnates  made  of  glass-epoxy  and  graphite-epoxy  com-  - 
posltes. 

Bader,  et  al.  [2]  attributed  the  layer-thickness  effect  on  transverse 
cracking  to  the  constraining  actions  from  the  adjacent  0*- layers.  The  pres¬ 
ence  of  the  0*-layers  not  only  arrests  the  crack,  but  also  restrains  it  from 
forming  in  the  first  place.  Depending  on  the  thickness  of  the  90*-layer,  the 
quantity  which  creates  transverse  cracks  is  associated  with  the  amount  of 
strain  energy  stored  in  the  90*-layer. 

The  exact  kinematics  in  the  transverse  cracking  formation  are  not  clear. 

A  closer  examination  of  the  cracking  process  indicates  that  the  crack  generally 
_ follows  the  fiber-matrix  interface,  when  viewed  in  the  plane.  For  example. 

Fig.  1.2  shows  an  x-ray  plane-view  of  a  [+25/902la  graphite-epoxy  laminate. 

It  is  seen  that  all  the  transverse  cracks  span  the  entire  width  of  the  test 
specimen.  Similarly,  a  photo-micrograph  in  the  thickness  view,  is  shown 
in  Fig.  1.3  for  the  same  laminate.  Here,  the  transverse  cracks  span  across 
the  thickness  of  the  90*-layer  and  terminate  at  the  90/25  interfaces.  These 
photographs  do  not  reveal  the  actual  growth  path  of  the  transverse  cracks. 
Rather,  it  tends  to  suggest  a  sudden  formation  of  a  large  plane  crack,  con¬ 
fined  by  the  thickness  of  the  90*-layer  and  the  width  of  the  specimen.  The 
dynamic  nature  in  the  crack  formation  requires  further  verification,  however. 


Examination  of  the  transverse  crack  surface  by  scanning  electronic 
microscope  reveals  a  generally  rough  failure  surface  scattered  with  tiny  epoxy 
debris  and  some  broken  fibers.  Clearly,  the  cracking  path  is  not  entirely  confined 
in  the  matrix  phase  of  the  composite. 

Let  us  now  turn  to  the  problem  of  free  edge  delamination  which  is  another 
frequently  observed  sub- laminate  fracture  mode  in  epoxy-based  laminates.  Pig. 

1.4  illustrates  schematically  a  free  edge  de lamination  crack  in  a  straight-edge 
specimen  under  uniaxial  tension.  It  shows  again  a  plane  crack  which  forms 
along  the  free  edge  of  the  laminate  and  propagates  inward  along  an  Interface 
of  two  adjacent  layers. 

The  driving  force  for  the  free  edge  delamlhation  crack  stems  from  the- 
Interlaminar  stresses  that  exist  near  the  free  edge  region  of  the  laminate. 

These  stresses  are  generally  singular  in  nature  with  high  stress  gradient.  As 
the  cracking  forms  and  propagates  inward,  the  singular  Interlaminar  stress  field 
moves  also  with  the  crack-tip.  The  process  is  apparently  stable  a a  is  evidenced 
by  numerous  experiments.  On  the  other  hand,  there  is  no  crack-arresting  mecha¬ 
nisms  in  the  process;  and  the  crack  can  therefore  propagate  indefinitely  until 
the  laminate  is  separated  through  the  cracking  interface.  At  this  point. 


Fig.  1.5  shows  ths  x-ray  plane-view  pictures  of  a  1+45/0/90]  graphite- 
epoxy  laalnate  undergoing  free-edge  de lamination  at  three  different  loading 
stagas,  where  a  stable  crack  growth  is  Indicated.  The  interface  In  which  the 


crack  is  supposed  to  be  contained  should  be  the  mid-plane  of  the  laminate. 
However,  a  side-view  micrograph  of  the  same  laminate  taken  at  the  loading  stage 
(b)  is  shown  In  Fig.  1.6.  It  is  seen  here  that  the  edge  crack  Is  actually 
formed  within  the  90*-layer,  but  it  zig-zags  along  the  length  of  the  specimen. 
Apparently,  the  crack  is  not  always  confined  in  any  one  given  layer  Interface. 


(»)  0>)  (c) 


Figure  1.5  Z-ray  Pictures  of  a  [+45/0/90],  Graphite-Epoxy 

Laalnate  under  Going  De lamination  in  the  Mid-Plane 


Micrograph  Shoving  Free  Edge  Delamination  in  the  90  layer 
of  a  Graphite-Epoxy  (+45/0/90)  Laminate 


Another  side-view  micrograph  of  this  laminate  taken  at  the  loading  stage 
(e)  is  shown  in  Fig.  1.6.  Of  particular  interest  is  the  formation  of  other 
dalaaination  cracks  in  the  ±45 '-interfaces  and  the  0/45  interfaces.  These  cracks 
are  generally  in  shear-dominant  modes.  Clearly*  the  growth  process  of  free  edge 
de lamination  is  much  more  coup lice tad  than  the  transverse  cracking  process  that 
has  been  discussed  earlier. 

Analytical  investigations  of  the  free  edge  dalaaination  problem  have  been 
mainly  efforts  for  calculating  the  boundary  layer  interlaminar  stresses.  Almost 
without  say  exception,  all  the  calculation  methods  were  based  on  the  effective 
property  assumption,  that  the  material  layers  in  the  laminate  are  Individually 
homogeneous  (but  anisotropic)  media.  Within  this  assumption,  material  and 
geometrical  discontinuities  exist  only  across  the  layer  Interfaces. 

Even  at  this  dimensional  level,  the  elasticity  boundary-value  problem  is 
quite  complicated;  only  numerical  solutions  for  the  free  edge  stress  field  have 
so  far  been  obtained  [3-6]*. 

The  knowledge  of  the  free  edge  stresses  has  provided  qualitative  correla¬ 
tions  with  some  delamination  experiments.  For  example,  in  e  series  of  strength 
tests,  Bjeletlch,  et  al.  reported  the  average  tensile  strength  of  six  families 
of  quasi- iso tropic  laminates  by  alternating  the  stacking  sequence  of  the  0*, 

90*,  and  ±45*  layers  [9].  An  edge  stress  analysis  identified  correctly  that 
the  laainate  of  [0/90/±45]#  developed  a  co^resslva  o^  stress  along  the  free 
edge.  Consequently,  the  laainate  suffered  no  de lamination  until  final  failure. 


*  Recently,  Wang  and  Choi  [7]  have  determined  the  elasticity  solution  for  the 
singularity  of  the  edge  stresses.  Other  efforts  in  computing  the  free  edge 
stresses  are  discussed  more  extensively  in  a  survey  paper  by  Sslmmon  [8]. 
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Figure  1.6  Micrograph  Showing  Delamination  Cracks  at  Several  Interfaces  at  Late 
Stage  of  Loading  of  a  Graphite-Epoxy  [+45/0/90]  Laalnate 


On  the  other  hand,  the  laminate  of  (+45/Q/90]g  construction  developed  e 
tensile  a  stress  along  the  free  edge  boundary*  which  caused  prematured 
delamination;  and,  the  growth  of  the  delamination  resulted  in  a  lowered 
laminate  tensile  strength. 

In  a  series  of  tests  using  graphite-epoxy  laminates,  lodini,  at  el. 

[10]  found  that  the  critical  stress  at  the  onset  of  edge  delamination  in  a 
[4-45  /0  /90  ]  ,  n  ■  1,2,3  family  varied  greatly  with  the  value  of  n.  The 
critical  stress  decreases  at  tha  rate  of  about  /n.  In  addition,  the  location 
and  the  mode  of  crack  also  varied  with  n.  Clearly,  the  actual  thickness  of 
the  material  layers  in  the  laminate  influences  the  mechanisms  of  delamination. 
Tet,  a  free  edge  stress  analysis  for  the  laminates  would  yield  identical  edge 
stress  field  for  all  values  of  n.  Thus,  e  stress  analysis  alone  cannot  explain 
the  initiation  and  growth  events  in  the  delamination  process. 

The  layer-thickness  effect  on  edge  delamination  is,  coincidentally, 
similar  in  nature  with  that  found  in  the  transverse  cracking  processes.  In 
both  cases  it  requires  e  rigorous  physical  end  mathematical  explanation. 

All  of  this  points  to  the  need  for  a  fracture  analysis  of  the  problems. 

In  particular,  the  critical  conditions  under  which  a  crack  initiates  and 
propagates  must  be  defined.  And  the  growth  behavior  as  e  function  of  loading 
also  needs  e  closer  examination. 

As  discussed  earlier,  transverse  cracks  and  delamination  are  sub-laminate 
cracks,  which  may  be  analysed  at  two  different  dimensional  levels.  At  the 
microscopic  level,  the  basic  assumption  is  that  in  each  material  layer  there 
is  a  distribution  of  micro-defects.  And,  a  similar  micro-defects  distribution 
may  also  exist  on  the  layer  interfaces.  Thus,  increasing  the  thickness  of  the 
material  layer  increases  the  probability  of  crack  initiation  and  growth. 
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This  reasoning  can  certainly  explain,  albeit  qualitatively,  the  observed 
layer- thickness  effect  on  transverse  cracking  f  2J  and  edge  delaninatlon  [10] . 
But  the  analysis  wist  include  the  interrelationships  bateau  the  stress  fields, 
the  defects,  defect  interactions  and  their  coalescence  aechanisas  at  the 
aicroscopic  level.  This  approach,  however,  la  practically  unattainable. 

The  alternative  is  to  perform  a  fracture  analysis  at  the  aacroscoplc 
level,  where  a  stress  analysis  based  on  layer-elasticity  nay  be  coupled  with 
the  classical  fracture  mechanics.  Within  this  context,  one  must  perfora  a 
number  of  fundamental  preparations  in  order  to  construct  a  physically  realistic 
fracture  model.  In  particular,  soma  conceptual  and  technical  questions  must 
be  resolved  first. 

As  in  any  analytical  model  construction,  idealisation  and  s Implication 
are  made  in  order  to  reduce  mathematical  complexity.  In  the  case  of  free- 
edge  delamination,  for  instance,  one  may  Idealise  the  crack  growth  as  a  self- 
similar,  or  co-planar  crack  contained  in  an  interface  between  two  layers. 

The  growth  is  one-dimensional  and  the  crack  front  is  represented  by  a  point. 
Similarly,  transverse  cracking  may  also  be  idealised  as  a  self -similar  crack 
along  the  matrix-fiber  interface  In  the  90*-layer.  In  either  case,  the  crack 
surface  is  assumed  flat  and  smooth.  In  reality,  of  course,  the  crack  surface 
is  ragged  and  warped  depending  on  the  exact  nature  of  the  crack,  even  whan 
viewed  at  the  macroscopic  level.  But,  the  idealisation  may  be  tolerated  In 
the  overall  modeling,  if  the  material  property  (e.g.,  the  critical  energy  re¬ 
lease  rate,  Ge  or  K^)  is  measured  as  a  bulk  quantity,  and  the  site  of  crack 
is  large.  This  kind  of  idealisation  is,  on  one  hand,  a  mathematical  necessity; 
and  on  the  other  hand,  it  wist  be  adjudicated  by  physical  experiment. 


The  above  consideration  brings  upon  another  conceptual  question*  which 

is  how  G  or  K  should  be  physically  determined.  Although  transverse  cracks 
c  c 

and  edge  delanination  nay  all  be  regarded  as  matrix-dominant  cracks*  their 
individual  resistance  against  crack  growth  when  aeaaured  aacroscopically  could 
differ  considerably.  Consider*  for  example*  the  two  different  crack  actions 
in  a  unidirectional  laminate  as  shown  in  Fig.  1.8.  Case  (a)  Illustrates  a 
mode -I  delanination  action  between  0*/0*  layers*  and  case  (b)  depicts  a  mode-1 
crack  action  between  90° /90*  layers.  The  two  crack  actions  may  produce  distinc¬ 
tive  crack  surfaces  at. the  microscopic  level*  and  thus  give  different  values 

for  G  which  is  measured  at  the  macroscopic  level, 
c 

In  a  series  of  teats,  conducted  recently  at  Texas  A  &  M  University  [11,12]* 
different  values  for  the  critical  energy  release  rate,  Gc*  were  found  for  a 

unidirectional  composite,  depending  on  how  the  crack  is  propagated.  In  particu¬ 
lar,  for  the  AS-3502  composite  system,  Gfi  measured  under  0°/0°  delamination 
action,  case  (a)  Fig.  1.8,  was  about  0. 8-0.9  lb/in;  while  G.,  under  90° /90* 
cracking  action,  case  (b).  Fig.  1.8,  was  found  to  be  1.2-1.35  lb /in.  Gc  for 
the  neat  resin  (without  reinforcement)  had  a  value  of  only  0.4  lb/in.  The 
differences  all  stem  from  the  Individual  characteristics  ot  the  crack  surfaces. 
This  clearly  raises  the  question  of  whether  G  can  be  regarded  as  a  general 

C 

material  property  for  the  class  of  the  so-called  matrix-dominated  cracks. 

The  difficulty  comes,  obviously,  from  the  fact  that  Gfi  is  measured  without 
regard  to  the  microscopic  details  of  the  cracking  surfaces. 

Closely  related  to  the  question  of  Ge>  there  is  also  some  concern  about 
a  possible  non-linear  yielding  that  occurs  in  the  vicinity  of  the  crack-tip, 
especially  under  shearing  crack  actions.  In  a  recent  experiment  by  Vandarklay 
[13],  it  wss  shown  thst  the  total  energy  release  rate  measured  under  mixed- 
mode  (G  )  creek  action  is  much  larger  than  under  a  purs  mods-I  action. 

Xell 


(b)  Tranavaraa  Cracking  Actlon(90°/90°) ,  Mode-I 


Figure  1.8  Schaaatlcs  of  Two  Different  Cracking  Action* 
In  Onidirectionel  Leainate 


This  suggests  that  energy  dissipation  In  the  mixad-node  cracking  may  contain 

some  amount  of  Inelastic  deformation  energy*  In  view  of  the  feet  that,  the 

epoxy  matrix  generally  exhibits  a  considerable  inelasticity  under  sheering 

deformation ,  the  measured  value,  of  G  under  nixed  mode  crack  action  may-  depend 

c 

on  the  amount  of  shearing'  action.  .  -  -  •  - 

Clearly,  all  these  questions  require  in-depth  investigation  before  a 
reliable  fracture  model  can  be  formulated. 

It  is  therefore  the  purpose  of  this  research  to  conduct  a  detailed  and 
comprehensive  examination  into  the  growth  mechanisms  of  transverse  cracking 
and  free  edge  delamination.  The  main  objective  is  to  formulate  a  predictive 
fracture  model  which  is  capable  of  describing  the  basic  behavior  in  the  two 
types  of  sub-laminate  crack  processes.  To  this  end,  the  macroscopic  approach 
based  on  layer-elasticity  is  taken,  and  the  physical  concepts  of  the  classical 
fracture  mechanics  will  be  employed.  In  order  to  resolve  some  of  the  funda¬ 
mental  questions  associated  with  this  approach,  an  extensive  experimental 
study  is  also  conducted,  with  results  correlated  with  the  fracture  analysis. 

Section  II  details  the  fracture  model  development.  A  review  of  the  basic 
methods  in  classical  fracture  mechanics  is  included.  And,  the  special  method 
of  the  crack-closure  integral  is  singled  out  for  adaptation  in  a  finite  element 
computational  scheme.  In  this  section,  the  crack  initiation,  growth  and 
growth  stability  criteria  for  transverse  cracking  and  free  edge  de lamination 
are  proposed. 

Section  III  contains  the  results  and  other  data  documented  in  the 
experimental  case  study.  It  includes  the  mecroscoplcally  measured  material 
properties,  crack  events  in  the  [0/90  1  ,  n  -  1,2, 3,4  family,  crack  events  in 
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the  [425/90  ],  n  ■  1/2,1, 2,3,4, 6, 8  family,  and  the  crack  events  in  the  [+45  1 

a  w  D  ® 


a  ■  2,3  £  tally.  In  particular,  in  tha  [+25/90n]g  £  sally,  both  eransvaraa 
cracking  and  edga  delaalnation  occur  either  independently  or  Interactively, 
depending  on  the  value  of  n.  Other  peculiar  aspects  in  the  crack  events  are 
staamarixed  and  discussed. 

Numerical  modeling  based  on  the  methods  developed  in  Section  II  are 
conducted  for  the  laminates  tested  In  Section  III.  Thus,  a  complete  cor¬ 
relation  between  the  analysis  and  experiment  is  presented  is  Section  IV. 

This  includes  the  modeling  of  the  transverse  cracking,  free  edge  del, ami nation 
and  their  interactions  effect. 

A  concluding  discussion  on  the  gains  as  well  as  shortcomings  in  the 
fracture  model  is  given  in  Section  V. 
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II.  DEVELOPMENT  OF  FRACTURE  MODELS 

• 

In  the  previous  section,  we  have  discussed  some  aspects  of  the  physical 
process  of  transverse  cracks  and  edge  da lamination  In  epoxy-based  laminates. 

It  is  found  that  the  fracture  events  cannot  be  explained  from  the  standpoint 
of  the  concept  of  strength  as  a  constant  material  property.  Actually,  there 
is  a  broad  range  of  failure  processes  for  which  the  concept  of  strength  la 
not  applicable.  This  fact  is  particularly  veil  known  for  brittle  failures. 

The  essence  of  the  classical  fracture  mechanics  is  not  to  apply  the  strength 
concept,  but  base  the  theory  on  the  actual  process  of  fracture.  It  Is  founded 
on  the  work  of  Griffith  [14],  published  originally  in  1920. 

2.1  The  Griffith’s  Energy  Release  Rate  Concept 

Griffith's  original  work  [14]  Is  concerned  with  a  thin,  brittle  plate 
which  Is  uniformly  stretched  in  one  direction  by  o.  Fig.  2.1.  The  plate  has 
a  through-crack  of  length  2a,  orientated  transverse  to  the  direction  of  the 
applied  stress  o.  The  length  of  the  crack  Is  assumed  small  compared  to  the 
dimension  of  the  plate;  a  «  L;  a  «  W.  Griffith  Introduced  the  surface 
energy  of  the  brittle  material  and  formulated  the  criterion  according  to  which 
the  existing  crack  begins  to  propagate  unstably.  He  postulated  that  an  In¬ 
crease  in  the  crack  size  causes  a  decrease  In  the  stored  strain  energy  near 
the  crack- tip;  and  the  loss  of  the  strain  energy  Is  converted  entirely  Into 
surface  energy.  The  latter  is  the  result  of  an  Increase  In  the  free  surface 
area  due  to  the  crack  extension. 

Let  U  and  UQ  be  the  total  strain  snergy  of  the  plate  with  and  without  the 
crack,  respectively.  And,  let  AU  ■  Dq  -  D  be  the  loss  of  strain  energy  of  the 
plate  due  to  the  presence  of  the  crack.  Then,  for  the  existing  crack  of  slse 
2a,  the  Griffith  criterion  states  that  the  crack  begins  to  propagate  when 
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r-  (AD)da  >  (4y)da 


(2.1) 


where  y  la  the  material  surface  energy  of  a  unit  of  free  surface. 

For  the  Griffith  problem  (plane  stress)*  AD  la  given  by 

AD  -  |  (2.2) 

Dpon  substitution  Into  (2.1)»  the  following  dependence  of  the  critical  stress  on 
crack  site  Is  obtained: 


According  to  Eq.  (2.3)*  the  crack  of  original  length  2a  remains  unchanged 

until  the  applied  stress  a  reaches  the  value  given  by  (2.3);  at  this  Instant 

a  dynamic  process  of  crack  growth  begins. 

The  Griffith  equation  (2.3)  brings  upon  two  Important  practical  questions. 

The  first  is  that  using  (2.3),  a _ ♦  •  is  a  +  0.  Of  course*  no  real  material 

cr 

can  sustain  an  Infinite  stress.  In  fact,  given  a  plate  having  no  crack  at  all* 
one  finds  only  a  finite  strength  for  the  plate.  This  contradic- *?.<*•  Is  cicum¬ 
vented  by  Introducing  the  concept  of  Inherent  material  flaws,  that  ars  assumed 

to  exist  naturally  within  the  plate.  The  real  physical  identity  of  the  flaws 
Is  lost  within  the  frame-work  of  the  "continuum  assumption"  from  which  stresses* 
strain  energy*  etc.  are  calculated.  In  reality*  the  flaws  exist*  perhaps*  at 
the  dimensional  level  of  the  grain-boundary,  or  the  molecular  structure  of  the 

material.  The  macroscopic  effect  of  these  mlcroflavs  is  then  represented  by 
assuming  a  macroscopic  crack  of  size  Sq*  from  which  the  no-crack  strength  of  the 
plate  can  be  calculated  from  (2.3).  Clearly*  the  existence  of  a^  at  the  macro¬ 
scopic  level  can  only  be  viewed  as  a  necessary  postulation. 
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Another  question  is  related  to  the  definition  of  f,  the  free  surface 
energy.  Actually,  it  represents  the  irreversible  work  required  to  create  a 
unit  free  surface  area.  For  brittle  material  such  as  glass ,  crack  growth  as 
defined  by  (2.3)  has  been  shown  to  be  quite  valid.  For  most  other  structural 
materials,  die  crack-tip  region  usually  exhibit  some  degree  of  ductile  deforma¬ 
tion.  Early  studies  by  Irwin  [15],  Orowan  [16]  and  others  led  to  the  concept 
of  the  so-called  "quasi-brlttle  fracture."  According  to  this  concept,  Eq.  (2.3) 
is  also  correct  for  most  materials  under  conditions  of  quasi-brlttle  fracture, 
if  the  quantity  y  is  replaced  by  the  irreversible  energy  dissipated  In  the 
surface  region  of  the  crack  per  unit  area  of  free  surface.  This  last  quantity 

depends  on  the  ductility  or  inelasticity  of  the  material  near  the  crack-tip; 
and  it  is  usually  found  to  be  orders  of  magnitude  larger  than  the  theoretically 
calculated  value  of  the  material's  free  surface  energy  y. 

The  foregoing  discussions  become  practically  more  important  when  we  try 
to  describe  crack- like  failures  in  composite  materials.  The  fundamental 
assumptions  in  the  Griffith  theory  must  be  propertly  re-lnterpretated  in  order 
for  it  to  apply  to  the  fracture  problems  of  composites. 

As  has  been  discussed  in  Section  I,  practical  stress  analysis  for  multi¬ 
layer  laminates  can  be  performed  only  at  the  level  of  layer-elasticity.  That 
is  the  individual  layer  is  idealized  as  a  homogeneous  medium  in  much  the  same 
manner  as  the  "continuum"  assumption  in  the  theory  of  elasticity.  Consequently, 
when  we  analyze  a  crack-like  failure  in  this  idealized  material,  the  same  con¬ 
cern  arises  in  identifying  properly  the  quantity  aQ  and  the  quantity  y  if  the 
Griffith  theory  is  to  be  employed.  As  before,  these  quantities  must  be  con¬ 
sidered  a  material  property  and  are  measured  at  the  dimensional  level  where 
the  stress  analysis  is  performed.  As  far  as  composite  materials  are  concerned. 


the  aiM  of  *0  can  ba  vary  larga,  depending  on  tha  slsa  and  tha  distribution 
of  tha  nicro flaws,  Similarly  tha  quantity  y  haa  to  ba  tha  averaged  value 
over  a  larga  crack  surface  area. 

In  this  sanaa*  tha  Griffith  aquation  (2.3)  can  ba  generalised  to  treat 
tha  types  of  crack  growth  problems  which  are  found  In  layered  coapoaltes. 

Now,  consider  tha  quantity  AU  in  (2.1).  Since  Oq  does  not  depend  on  ar 
the  Griffith  criterion  may  be  rewritten  as 


30  3S. 

3a  3a 


(2.4) 


where  S  is  the  total  surface  energy  of  the  solid  In  the  context  of  the  afore¬ 
mentioned  "quasi-brittle"  fracture. 

The  quantity  3U/3a  depends  on  the  stress  field  near  tha  crack-tip,  which 

in  turn  depends  on  the  geometry  of  the  crack.  It  represents  the  driving  force 

in  propagating  the  crack  and  is  commonly  referred  to  as  the  available  strain 

release  rate,  G(a),  a  function  of  crack  size  a.  The  quantity  3S/3a  depends 

-  on  the  imcroatructure  of  the  crack  surface  and  it  represents  the  resistance 

of  the  material  against  the  crack  propagation.  For  material  of  uniform 

property,  3S/3a  is  independent  of  a.  It  is  commonly  called  the  critical 

energy  release  rate  G  ,  which  is  identified  as  2y. 

c 

Thus,  the  Griffith  criterion  (2,4)  becomes, 

G(a)  £  Gc  (2.5) 

Accordingly,  the  development  of  the  classical  fracture  mechanics  rests 
upon  the  calculation  of  G(a)  analytically,  and  the  measurement  of  Gfi  physically. 
In  view  of  the  mathematical  complexity  in  the  singular  stress  field  near  the 

crack-tip  and  the  practical  uncertainties  involved  in  the  material  behavior  at 
failure,  both  tasks  have  not  been  a  simple  matter. 
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2.2  Irwin’s  Crack-Closure  Integral 

For  plane  problems  in  elastic  hoaogeneous  solids,  the  strain  and  dis¬ 
placement  variables  are  Independent  of  one  of  the  Cartesian  coordinates ,  say 
z.  In  this  general  case,  the  singular  stress  field  near  the  crack  tip  can  be 
represented  by  analytical  functions  in  the  theory  of  complex  variables  [17]. 
Generally,  stress  field  solutions  are  obtained  for  three  particular  nodes  of 
crack  action.  These  three  modes,  illustrated  in  Fig.  2.2,  are  known  respec¬ 
tively  as  mode  I  or  opening  mode,  mode  II  or  sliding  mode  and  mode  III  or 
tearing  mode.  For  isotropic  material,  the  distribution  of  stresses  near  the 
tip  of  an  arbitrary  brittle  crack  undergoing  pure  mode  I  extension  is  given 
by  [18]. 
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2.2 


And  the  crack-tip  stress  field  of  a  node  II  crack  is 
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And  the  crack-tip  stress  field  of  a  mode  III  crack  is 
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The  real  parameters  K^f  and  Kjjj  are  known  as  the  stress  intensity 
factors*  which  depend  generally  on  the  shape  of  the  body,  loading*  location 
and  size  of  the  crack. 

Knowing  the  crack-tip  stress  field  and  the  functional  dependence  of  the 

stress  Intensity  factors  on  the  crack  size  a,  the  corresponding  value  of  the 

available  energy  release  rate  G(a)  can  be  calculated. 

,2 

G(a)  ■  Kj  (plane  strain)  (2.9) 

or 


G(a)  -  ^  Kj2  (plane  stress) 
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(2.10) 
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In  the  general  case  of  mixed-mode  creek  action,  the  energy  release  rate 


Is  given  by 


G(a)  - 


1-v  _  2 


2V  .  1+v  2 


+  *ii  > +  «  *in 


(2.11) 


Thus,  by  applying  the  Griffith  criterion  (2.5) ,  we  relate  G  as  follows: 

c 


sc-T^«ic2  +  *iic2>+TL*iii«2 


In  particular,  for  pure  mode  I  crack  extension, 

2 

GIc  -  Kjc2  (plane  strain) . 


(2.12) 


And  similarly. 


lie  E 


,  -  3±v 

*IIIc  E 


(2.13) 


(2.14) 


(2.15) 


Fran  the  original  statement  of  Griffith,  we  can  identify  only 
G^  ■  Gjc  “  2y.  For  mixed-mode  cracks,  experience  has  shown  that  the  value  of 
G^  depends  on  the  ratios  of  K^/IL^  and  In  fact.  In  most  brittle 


eats  rials,  G  is  much  less  titan  G  fovnd  under  mixed-mode  actions. 

c 

Thus,  traditionally  some  seedL-eaplrlcal  crack  growth  criteria  are  used. 


depending  on  the  particular  situation.  These  criteria  are  In  the  general  form 


f<*ic-  *nc’  *IUC>  -  °- 


(2.16) 


Mote  that  In  a  criterion  such  as  (2.16),  the  quantities  K  *s  and  G  'e 

c  c 

are  Interchangeable  through  the  one-to-one  relatione  (2.13)  to  (2.15). 

Relations  between  K  and  G  for  orthotropic  homogeneous  media  with  the 
crack  orientated  along  one  of  the  aajor  axes  can  also  be  obtained  [19].  But, 
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for  layered  aadla  such  as  composite  laminates,  an  explicit  relationship 
between  K  and  G  Is  usually  not  available.  Moreover,  in  a  generally  aniso¬ 
tropic  solid,  the  mode  I,  mode  II  and  mode  III  crack  actions  are  mutually 
coupled. 

Sven  with  Isotropic  msdla,  tha  evaluation  of  the  stress  intensity  factors 
may  become  difficult  if  not  untractable.  Irwin  [20]  showed  that  the  elastic 
strain  energy  released  during  an  Incremental  crack  extension  can  be  equated 
to  the  work  done  in  closing  the  Incremental  crack.  In  what  Is  known  as  the 

crack-closure  representation,  the  method  provides  for  a  direct  determination 
of  G's  from  the  crack-tip  stresses  and  displacements. 

Let  Aa  represent  an  Infinitesimal  crack  extension  from  the  initial  sice 
of  a  as  shown  in  Figure  2.3.  Let  Au  be  the  relative  displacements  between 
the  mating  crack  surfaces  along  Aa  when  the  crack  is  open  as  In  Figure  2.3b. 
Let  o  be  the  surface  stress  distribution  along  Aa  when  the  crack  is  closed  as 
In  Figure  2.3c.  Then  the  work  done  to  close  the  crack  extension  Is  given  by 

rAa 

AW  -  1/2  a  •  Au  da  (2.17) 

J0 

Thus,  by  means  of  Irwin's  crack  closure  equivalence,  the  available  energy 
release  rate  G,  for  a  crack  of  size  a,  Is  expressed  by  the  integral, 

“<*>  *  24l  f0‘  5  •  41  **  (2-18> 

Substituting  the  components  of  the  surface  stresses  3  and  the  relative  dis¬ 
placements  Au  into  equation  (2.18)  yields  G  in  component  form  as 
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where  G^,  Gjj,  G^^  denote  the  node  I,  node  II  and  node  III  components  re¬ 
spectively;  and  the  total  energy  release  rate  under  generally  mixed-mode  crack 
action  Is, 

G(a)  -  Gx  +  Gu  +  GIH.  (2.20) 

The  crack-closure  representation  Is  particularly  convenient  for  adapta¬ 
tion  In  numerical  computation.  Other  methods  such  as  the  well  known  J-integral 
technique  [21,22]  are  also  used  to  evaluate  G(a).  These,  however,  will  not 
be  detailed  here,  as  they  are  not  pertinent  to  our  subsequent  development. 

2.3  Finite  Element  Representation 

A  numerical  technique  to  calculate  the  strain  energy  release  rate  G(a) 
has  been  presented  by  Ryblcki  and  Kanninen  [23] .  Their  approach  Involves  a 
finite  element  solution  of  Irwin's  crack  closure  integral  given  in  (2.18). 

In  the  finite  element  representation,  the  continuous  stress  and  displacement 
fields  of  the  solid  are  approximated  by  the  nodal  forces  and  displacements, 
respectively.  Fig.  2.4  Illustrates  the  finite  element  representation  of  a 
crack  tip  region.-  Here,  a  crack  of  length  a  is  shown  with  the  crack  tip  at 

node  c.  The  finite  element  solution  determines  the  displacement  components 

(u,v,w)  of  the  crack  tip  node  c.  An  Incremental  crack  extension  Aa  Is  lntro- 
c 

ducad  by  replacing  the  crack  tip  node  c  with  two  separate  nodes  f  and  g  as  shown. 


29 


*7 


With  this  new  crack  geometry  taken  into  account,  the  finite  element  solution 

for  the  nodal  displacements  (u,v,w)_  and  (u,v»w)  are  found  for  nodes  f  and 

*  © 

g  respectively.  The  crack  extension  is  then  closed  by  applying  equal  and 
opposite  forces  at  nodes  f  and  g  such  that  their  common  displaceaents  aatch 
the  displaceaents  found  earlier  for  node  c.  Of  course,  these  forces  are 

actually  the  Internal  nodal  forces  which  existed  at  node  c  before  it  is  opened. 
The  work  required  to  close  the  crack  extension  is  approximated  by 

AW~[Fx(uf-ug)  +  Fy(vf-vg)  +  Fz(wf-wg)]/2  (2.21) 

where  F  ,F  ,F  are  the  components  of  the  nodal  forces  required  to  close  nodes 
x  y  z 

f  and  g  together.  Thus,  the  energy  release  rates  for  the  three  crack  extension 
modes  are  approximated  by 

Gi  ~  Fy(vf~vg)/2Aa 

Gt_  ^  F  (u.-u  )/2Aa  (2.22) 

11  x  r  g 

GIU  ~  VVV/2Aa 

The  method  presented  here  does  not  require  that  the  stresses  be  calcu¬ 
lated  because  the  stress  and  strain  fields  are  approximated  by  the  nodal  forces 
and  displaceaents  in  the  finite  element  solution.  Rybicki  and  Kanninen  (23 2 
applied  conventional  constant  strain  finite  elements  to  the  solution  of  three 
fracture  problems:  the  double  cantilever  beam  specimens;  a  finite  strip  con¬ 
taining  a  central  crack;  and  a  bolt  fastened  double  lap  joint  containing  radial 
cracks.  They  reported  that  a  relatively  coarse  grid  in  the  crack  tip  region 
was  sufficient  for  a  good  comparison  between  their  example  problems  and  refer¬ 


ence  solutions. 


Nevertheless,  Che  adequacy  of  this  kind  of  approximation  remains  a  sub¬ 
ject  of  great  concern  because  the  mathematical  singular  nature  of  the  crack- 
tip  is  disregarded  in  the  approximation.  But,  even  within  the  exact  definition 
of  6(a)  such  as  by  Irwin's  integral  (2.18),  the  quantity  3  •  4u  at  the  crack- 
tip  must  be  finite  (cr  is  singular,  Au  is  zero).  Thus,  the  approximation  of 
6(a)  in  (2.21)  is  logically  acceptable.  The  only  question  is  the  degree  of 
accuracy.  The  latter  is  ultimately  related  to  the  manner  under  which  the 

material  properties  are  determined. 

In  this  report,  the  numerical  technique  described  above  is  applied  to 

calculate  the  strain  energy  release  rate  as  a  function  of  crack  length  for  the 
transverse  cracking  and  delamination  fracture  processes  in  composite  laminates . 
The  formulation  of  the  actual  finite  elements  appropriate  in  these  calculations 
are  presented  in  Section  IV.  The  detiala  of  the  fracture  initiation  and  growth 
criteria,  however,  will  be  presented  next. 

2.4  Initiation  and  6rowth  Criteria 

In  this  study,  we  investigate  two  types  of  sub-laminate  cracks.  Namely, 
the  90°-cracking  and  the  free  edge  delamination.  In  order  to  set  up  a  crack 
initiation  and  growth  criterion  for  each  of  the  two  types  of  cracking,  we 
assume  that  micro-flaws  exist  along  the  fiber-matrix  interface,  and  the  ply- 
to-ply  Interfaces.  These  flaws  are  generally  of  a  size  in  the  order  of  the 
fiber  diameter,  whose  crack  propagation  mechanisms  is  in  the  realm  of  micro- 
mechanics.  Here,  we  postulate  only  that  these  micro-flaws  propagate  and 
coalesce  into  a  macroscopic  proportion  under  a  certain  critical  far  field  load¬ 
ing  condition;  the  initiation  process  is  dynamical  at  the  instance  of  forma¬ 
tion.  The  size  of  the  crack  at  the  macro-scale  depends  on  the  micro-structure 
of  the  material  (perhaps  also  on  a  host  of  other  material  processing  factors, 
such  as  the  curing  process,  post-cure  handling,  etc).  This  initial  macrocrack 
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is  called  here  Che  initial  "flaw"  whose  size  Is  denoted  by  a^.  The  value  of 
*q  is  generally  unknown  and  can  be  random  in  nature.  It  is  regarded  here,  how¬ 
ever,  as  a  known  macroscopic  material  property  along  with  other  properties 

such  as  the  stiffness  constant,  G  ,  etc.  of  the  basic  material. 

c 

As  has  been  discussed  previously,  both  the  transverse  cracking  and  free 
edge  delamination  processes  may  be  modeled  as  one-dimensional  self-similar 
crack  propagation.  Thus,  the  energy  release  rate  G(a)  in  each  case  can  be 
calculated  by  the  numerical  finite  element  technique. 

The  Griffith  criterion  of  (2.5),  when  applied  to  predict  the  initiation 
of  the  sub-laminate  cracks,  is  replaced  by  the  general  statement. 


S(.„>  >  G, 


•0-  -  c 

And,  the  stability  of  growth  following  initiation  is  governed  by  the 


(2.23) 


conditions , 


G(a  +  Aa)  <  G£ 

stable  growth 

(2.24) 

G(a  +  Aa)  >  Gc 

unstable  growth 

(2.25) 

G(a  +  Aa)  G£ 

neutral  growth 

(2.26) 

where  a  (>  aQ)  is  the  size  of  an  existing  crack  for  which  the  growth  stability 
is  to  be  determined.  We  must  note  here  that  the  definition  of  aQ  is  in  accordance 
with  the  original  Griffith  problem.  It  is  a  property  of  the  basic  material  system. 

(a)  Transverse  Cracking  Model 

The  crack  growth  geometry  for  a  single  transverse  crack  is  Illustrated  in 
Fig.  2.5.  Under  the  idealized  cracking  geometry,  as  shown  in  Fig.  2.5(a), 
assume  that  the  thickness  of  the  90s-layer  2b  >  2sq.  The  crack  growth  path 
from  the  initial  size  of  2a^  is  assumed  to  initiate  in  the  middle  of  the 
90°-layer.  The  growth  geometry  is  then  represented  by  the  finite  element  model 
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as  shown  in  Pig.  2.5(b).  Thus,  the  solution  domain  is  the  2-dlmenslonal 
x-z  plane,  resulting  in  a  generalized  plane  strain  deformation  field. 

Fig.  2.6  shows  a  typical  finite  element  mesh  for  the  growth  of  a  trans¬ 
verse  crack  in  the  upper  half  of  the  laminate.  If  the  laminate  is  subjected 
to  a  uniaxial  tension  represented  by  a  far  field  uniform  displacement  5, 
the  energy  release  rate  function  G(a)  may  be  calculated  and  expressed  by 


G(a)  -  [Ce(a)Sx2]t,  (2.27) 

where  t  is  a  characteristic  length  (e.g.,  thickness  of  a  ply)  and  is 

the  far  field  laminate  strain  induced  by  6.  The  coefficient  function  C  (a) 

e 

is  calculated  by  imposing  e^  ”  1.  Note  that  Ce(a)  is  Independent  of  the 
loading  6. 

Fig.  2.7  shows  the  general  behavior  of  (a)  for  the  crack  in  the 
laminate  whose  90°-layer  thickness  is  larger  than  2aQ.  Note  that  the  crack 
size  2a  is  limited  by  the  thickness  of  the  90°-layer,  2b.  The  strain  energy 
release  rate  increases  initially  with  the  Increase  of  cracking  size;  as  the 
crack  approaches  the  ply  interface,  the  energy  release  rate  begins  to  decrease 
the  rate  of  decrease  depends  on  the  relative  rigidity  of  the  outside  layer. 
Generally,  the  maximum  of  G(a)  is  located  at  a  <  b. 

SL 

Similarly,  if  the  laminate  is  subjected  to  a  uniform  temperature  drop 
of  At,  and  if  tensile  stress  in  the  90*-layer  is  also  induced,  then  the  calcu¬ 
lated  energy  release  rate  function  for  the  transverse  crack  is  expressed  by 

G(a)  -  {CT(a)AT2]t.  (2.28) 
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Figure  2.6  Typical  Finite  Element  Mesh  for  Transverse  Cracks. 


Typical  behavior  of  (a)  for  a  transverse  crack  Is  similar  to  that  of 

(*)  «s  shown  In  Fig.  2.7.  And,  C^(a)  has  a  maximum  value  also  at  the  same 

a  <  b.  The  entire  C_(a)  curve  Is  generated  numerically  by  Imposing  AT  -  -  1. 
m  i 

In  practice,  all  epoxy-based  laminates  are  subjected  to  some  thermal 
residual  loading  when  cured.  The  effects  of  the  residual  stresses  on  trans¬ 
verse  cracking  can  be  evaluated  by  simply  assuming  a  uniform  temperature  change. 

AT  -  T  -  Tq  (2.29) 

where  Tq  is  usually  the  curing  temperature,  or  the  temperature  at  which  the 
laminate  Is  free  of  residual  stresses. 

Thus,  when  the  laminate  Is  loaded  by  uniaxial  tension,  denoted  by  the 
far-field  strain  e^,  a  combined  loading  condition  must  be  considered.  In 
this  case,  the  total  energy  release  rate  6(a)  is  expressed  In  the  form 

6(a)  -  [Ce(a)Sx2  +  CeT(a)Sx  AT  +  CT(a)AT2]t  (2.30) 

where 

CeT(a)  -  2(Ce(a)-CT(a)]1/2  (2.31) 

The  growth  of  transverse  cracking  is  primarily  In  mode-I,  or  the  open¬ 
ing  mode.  Accordingly,  onset  of  a  transverse  crack  in  the  90°-layer  whose 
90°-layer  thickness  Is  larger  than  2Sq  Is  defined  when 

G(a0)  -  GIc,  if  b  >  aQ  (2.32) 

With  the  definition  (2.32)  for  crack  Initiation,  Eq.  (2.30)  becomes  a 

quadratic  equation  in  ex»  if  AT  is  given.  The  positive  root  of  ex  defines 

the  critical  applied  laminate  stress  at  onset  of  the  crack, 

(o  )  -  I  (e  )  (2.33) 

x  cr  x  x'cr 

where  §x  Is  the  laminate  stiffness  In  the  loading  direction. 
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Figure  2.8  Criterion  for  Onset  of  a  Transverse  Crack  In  90° -Layer 
with  Thickness  b  >  a.. 


Fig.  2.8  illustrates  the  criterion  for  the  onset  of  transverse  crack 
graphically.  It  is  seen  that  G(sq)  «t  tha  applied  laminate  stress  given  by 

(2.33)  has  a  value  equal  to  G„  .  Once  the  crack  la  f oread,  the  available 

lC 

G(a)  >  GIc  for  all  a  >  a^.  Bancs,  the  growth  is  unstable,  or  is  dynamical. 

Tha  crack  is,  of  course,  eventually  arrested  by  the  outside  constraining  layers. 
During  this  unstable  crack  growth,  excess  energy  is  released,  represented  by 
the  shaded  area  in  Fig.  2.8.  The  effect  of  the  excess  energy  on  the  laminate 
will  be  discussed  later  in  this  report. 


In  cases  where  the  thickness  of  the  90*-layar  is  small,  such  as  the 


[0/90/0]  laminate,  then,  b  ■  one-half  the  90° -ply  thickness.  Since  a^  is  a 
property  of  the  basic  material  determined  within  the  premise  of  original 
Griffith  problem,  while  b  can  be  made  arbitrarily  small,  it  is  possible  that 
b  <  Sq.  In  this  case,  the  avllable  energy  release  rate  is  limited  by  b.  Fig. 

2.9.  In  fact,  the  maximum  available  G  ia  at  a  -  a  <  b.  Hence,  the  onset 

a 

of  transverse  crack  is  defined  by 


6<V  "  Glc’  i£  b  <  V 


(2.34) 


Eq,  (2.34)  determines  essentially  the  smallest  possible  load  for  a  transverse 


Figure  2.9  Criterion  for  Onset  of  a  Transverse  Crack 
in  90*-layer  with  Thickness  b  <  an. 
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It  is  seen  that  the  ctack  ia>  started  at  the  size  of  a  *  a^;  and  the  growth 
Is  stable  because  the  available  G(a)  <  G£c  as  a  grows  larger  than  a^.  A 
slightly  higher  load  is  then  required  to  propagate  the  crack  to  the  interface. 

The  foregoing  discussion  is  concerned  with  the  critical  condition  under 
which  a  transverse  crack  is  formed.  If  the  macroscopic  properties  of  the  lami¬ 
nate  are  ideally  uniform  throughout,  then  theoretically,  multiple  transverse 
cracks  in  the  90a-layer  will  appear  simultaneously  at  a  regular  spacing.  The 
latter  is  sometimes  referred  to  as  the  characteristic  spacing  [24]. 

The  characteristic  spacing  of  multiple  transverse  cracks  has  been  explained 

by  the  so-called  shear-lag  [25].  It  is  simply  that  an  interlaminar  shear  stress 

txx  (refer  to  coordinates  in  Fig.  2.5)  is  developed  along  the  interface  where  a 

transverse  crack  is  terminated.  This  shear  stress  is  highly  concentrated  near 

the  root  of  the  crack  and  it  decays  rapidly  at  a  distance  away  from  the  crack 

root.  Fig.  2.10  shows,  for  example,  the  shear  stress  t  acting  along  the 

zx 


Figure  2.10  Interlaminar  Shear  Stress  on  -25/90  Interface. 

x  is  measured  from  the  root  of  the  crack. 
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-23/90  Interface  of  a  graphite-epoxy  [+25/90n]fl,  n  •  1*2,3. 

Fro*  Fig.  2.10,  it  is  seen  that  the  shear  stress  is  singular  at  the  crack 
root*  and  it  diminishes  completely  at  a  distance  of  about  4b  away  fro*  the 
crack  root.  As  the  interlaminar  shear  stress  diminishes,  the  tensile  stress  a ■ 
in  the  90*-layer  increases;  it  reaches  the  far-fleld  value  outside  the  shear- 
lag  zone.  A  second  transverse  crack  can  for*  only  outside  this  zone.  Hence, 
the  so-called  characteristic  spacing  between  two  neighboring  cracks  should  be 
equal  or  greater  than  the  size  of  the  shear-lag  zone.  In  the  case  of  the  ex- 
antple  Illustrated  in  Fig.  2.10,  the  characteristic  spacing  S  is  4  times  the 
thickness  of  the  90*-layer  (S  -  8b).  Or,  the  characteristic  crack  density 
(cracks  per  unit  length)  is  N  -  1/S. 

Indeed,  if  we  place  the  second  crack  outside  the  shear-lag  zone,  see 
Fig.  2.6,  and  generate  the  energy  release  rate  curve  G(a),  an  identical  curve 
will  be  obtained  as  that  for  the  first  crack.  On  the  other  hand,  if  we  place 
the  second  crack  inside  the  shear- lag  zone,  the  energy  release  rate  curve  would 
be  smaller  than  that  for  the  first  crack  [26]. 

In  reality,  the  Idealized  oultiple  cracking  process  does  not  occur. 

Rather,  if  the  first  crack  foras  at  soae  critical  far-field  load,  o  .  then 

cr 

the  second  crack  would  form  at  a  different  location  (not  necessarily  at  the 
characteristic  spacing)  under  a  slightly  higher  load,  and  so  forth.  Hence* 
the  density  of  the  cracks  is  generally  rising  with  the  applied  load,  such  as 
illustrated  in  Fig.  2.11.  Horeover,  no  characteristic  crack  density  can  be 
defined.  And,  often,  the  event  of  multiple  crack  formation  la  interrupted  by 
other  sub-laainate  cracking  event  at  higher  load,  such  as  edge  da lamination* 
or  other  interface  crackings. 

Assuring  the  formation  of  each  transverse  crack  is  the  result  of  an 


Initial  macro-flaw  in  the  90 '-layer,  it  may  thus  be  appropriate  to  postulate 
that  a  random  distribution  of  the  initial  macro-flaw  exists  in  the  90°-layer. 
But,  this  aspect  pertains  to  the  statistical  nature  of  material  propexty,  not 
the  mechanism  of  crack  formation.  Therefore,  we  shall  defer  the  subject  to  a 
later  study. 


Figure  2.11  Transverse  Crack  Density  vs.  Applied  Load 

N  is  the  Theoretical  Characteristic  Density 


(b)  Free  Edge  Delamination  Model 

The  cracking  geometry  of  free  edge  delamination  in  a  symmetrical  laminate 
is  illustrated  in  Figs.  2.12  and  2.13.  Fig.  2.12  shows  a  mid-plane  delamlna- 
tlon  and  the  corresponding  finite  element  modeling  scheme;  and  Fig.  2.13 
shows  a  delamination  along  an  interface  other  than  the  mid-plane  of  the  lami¬ 
nate.  Owing  to  the  symmetry  of  the  laminate,  the  cracking  action  in  the  former 
case  is  essentially  of  mode  I,  while  in  the  latter  case  the  crack  action  is 
generally  mixed,  Including  modes  I,  II  and  III. 


la  both  cases,  the  crack  Is  assumed  to  form  along  the  prescribed  inter¬ 
face,  and  to  propagate  In  self-similar  manner  toward  the  interior  of  the 
laminate.  Hence,  the  solution  domain  Is  the  2-dlmenslonal  y-z  plane,  being  a 
generalized  plane  strain  deformation  field. 

Fig.  2.14  shows  a  typical  finite  element  mesh  for  the  growth  model  of  an 
edge  crack  along  a  given  layer  Interface.  It  is  assumed  that  an  initial  macro¬ 
flaw  of  size  ag  exists  near  the  free  edge  and  it  can  be  propagated  once  the 
applied  far-field  laminate  stress  o^  reaches  a  certain  critical  value. 

In  order  to  determine  this  critical  condition,  the  energy  release  rate 
curve  G(a)  is  first  generated  numerically  by  the  finite  element  procedure. 

Since  a  laminate  has  more  than  one  interface,  it  is  not  generally  possible  to 
pin  point  exactly  which  interface  is  going  to  delaminate  under  the  applied  load. 
Normally,  a  free  edge  interlaminar  stress  analysis  is  conducted  first,  before 
any  edge  delamination  calculation.  The  details  of  the  interlaminar  stress  field 
can  provide  some  indication  as  to  on  which  one  of  the  interfaces  an  edge  crack 
is  likely  to  occur.  Then,  a  calculation  of  G(a)  for  crack  propagation  along 
the  selected  interface  is  performed  next.  In  case  there  are  more  than  one 

possible  interfaces  which  are  likely  to  delaminate,  then  the  G(a)  curves  for 

\ 

each  of  these  possible  cracks  have  to  be  generated. 

As  has  been  discussed  earlier,  mid-plane  delamination  is  primarily  of 
mode  I,  or  opening  mode,  in  a  symmetric  laminate;  while  an  off -mid-plane 
delamination  is  generally  a  mixed-mode  crack.  Between  the  two  different  modes 
of  delamination,  the  crack  growth  criteria  may  also  be  different.  Depending 
on  the  microscopic  details  of  the  cracking  surfaces,  the  critical  energy  release 
rates  for  mid-plane  cracking  and  for  off-mid-plane  cracking  may  actually  be 
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Figure  2.14  Typical  Finite  Element 


distinct.  And  this  has  been  shown  experimentally  for  some  unidirectional 
graphite-epoxy  [13].  Hence,  a  trial  calculation  is  usually  required  to 
determine  just  which  interface  is  the  most  energetically  possible  for  an  edge 
delamination. 

Fig.  2.15  shows  a  typical  6(a)  curve  for  a  delamination  growth  along 

either  the  mid-plane,  or  an  interface  other  than  the  mid-plane.  The  general 

behavior  of  G(a)  is  that  it  will  Increase  rapidly  to  reach  a  maximum  value 

at  a  »  a  .  After  this  point,  the  value  of  G  becomes  essentially  a  constant, 
m 

independent  of  the  crack  size  a  (assuming  the  width  of  the  laminate  Is  large 


a 


Figure  2.15  Typical  G(a)  for  Edge  Delamination. 
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compared  to  the  crack  size  a).  The  value  of  a_  at  which  G  attains  Its  maximum, 

m 

depends  on  a  number  of  geometrical  factors.  The  most  notable  factor  is  the 
layer  thickness  of  the  two  adjacent  layers  on  whose  interface  the  delamination 
propagates.  For  example,  if  G(a)  represents  the  total  energy  release  rate 
curve  for  a  mid-plane  delamlnatlon  for  the  laminates  tj^5n/0n/90n]g,  a  *  1*2,3, 
then  the  maximum  of  G  would  occur  at  about  a  -  nt,  t  being  the  thickness  of 

ft 

one  ply.  Since  G  Increases  with  a  before  a  -*•  a  ,  the  available  »««<*«  G  In 

ni 

the  case  of  n  •  1,  for  Instance,  is  much  smaller  than  the  available  maximum  G 


in  the  case  of  n  ■  2, 

This  thickness  dependence  nature  of  G  on  layer  thickness  is  not 

max 

exactly  the  same  as  the  previously  discussed  transverse  cracking  problems.  But, 
the  reason  for  the  exlstance  of  a  maximum  is  also  due  to  the  constraining  effect 


The  physical  reason  for  the  difference  between  (2.36)  and  (2.37)  Is  the  seas 
as  In  the  transverse  cracking  problems  discussed  earlier. 

If  the  computed  6(a)  Is  for  an  off-mid-plane  de lamination  growth,  then 
G(a)  so  computed  is  consisted  of  three  modal  parts,  as  Is  expressed  generally 
in  Eq.  (2.20).  A  general  criterion  for  the  onset  of  such  a  mixed-mode  crack 
can  thus  be  similarly  given  as 


0(V  '  “(I.II.IIDc’  l£  *0  *  *» 
“<*»>  -  “(X.IX.IIDc’  1£  *0  ”  *. 


(2.38) 

(2.39) 


In  actual  problems,  experiment  is  required  to  decide  whether  or  not 
GJc  Is  equal  to  G^  or  whether  or  not  G^  ^  depends  on  the 

ratios  of  Gj/Gjj  and  Gj/Gjjj,  «tc. 

Whichever  final  growth  criterion  .Is  used,  the  growth  stability  of 
delamination  is  theoretically  stable  if  a_  <  a:  and  it  Is  neutral  if 

u  Xu 

«0  >  aa.  In  reality,  however,  the  growth  is  generally  stable.  This  is  because 
(l)  the  Interface  flaw  size  a^  is  not  uniform  along  the  specimen;  the  weakest 
location  always  starts  delamination  first;  and  (2)  the  actual  width  of  the 
specimen  is  finite,  say  1  inch  wide;  as  the  delamination  growth  becomes  larger 
relative  to  the  width,  the  structural  stiffness  of  the  specimen  is  reduced; 
and  so  is  the  available  energy  release  rate  G  under  the  same  load. 

This  last  aspect  on  growth  stability  will  be  discussed  again  in  Sections 
111  and  IV  where  we  shall  examine  the  relevant  experimental  results  as  they 
are  correlated  with  their  computed,  counterparts. 
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III.  EXPERIMENTAL  STUDY 


•  3.1  Scone  of  Experiment 

The  main  objective  o£  the  experiment  is  to  conduct  an  extensive  case- 
study,  which  will  complement  the  development  of  the  analytical  fracture  models. 
As  has  been  discussed  In  the  previous  section,  there  are  several  physical  and 
conceptual  questions  that  must  be  resolved  by  experiment  if  the  methods  of 
classical  fracture  mechanics  are  to  be  used  in  the  fracture  analysis  of  sub- 
laminate  failures.  Thus,  the  experimental  study  will  provide  not  only  physical 
insights  in  model  development,  but  also  guidance  for  an  effective  numerical 
computation.  The  experiment  identifies  the  important  material  and  geometrical 
parameters  which  influences  the  various  fracture  processes. 

Accordingly,  the  experimental  study  consisted  of  the  following  major 
tasks: 

(*)  Basic  Material  Characterization.  Since  all  the  laminates  tested 
were  made  by  laminating  unidirectional  plies,  the  basic  ply  properties  were 
characterized  within  the  frame  work  of  "ply  elasticity,"  or  the  macroscopic 
"effective"  moduli  approach.  Both  the  elastic  compliances  and  the  strength 
properties  were  determined  by  experiments. 

In  addition  to  the  unidirectional  ply  properties,  same  basic  laminates 
were  also  tested.  Among  them  included  [Og] ,  [90g],  [+25]2g»  and 

[+45g]s.  These  were  tested  for  their  laminate  tensile  modulus  (Eg)  and  their 
laminate  tensile  strength  (ou)  only.  No  attempt  was  made  to  examine  their 
failure  mechanisms. 

The  proper  values  for  the  residual-stress-frea  temperature,  Tq,  has  also 
been  investigated  experimentally. 

i 

i 
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These  basic  material  data  ara  required  as  input  information  in  the  numer¬ 
ical  simulation  of  fracture  growth. 

(b)  Transverse  Cracks  in  [0/90^0]  -  Family t  This  series  of  tests  in¬ 
cluded  laminates  in  the  form  [0/90q/0],  n  -  1, 2,3,4  and  [O^O^/O^] ,  n  "  1»2, 
Because  of  the  relatively  small  edge  normal  stress  oc,  these  leminate  do  not 
suffer  edge  delamination  under  uniaxial  tension.  Rather,  only  transverse 
cracks  in  the  90°-layer  can  be  induced.  Hence,  this  series  of  testa  provides 
an  ideal  case  study  for  the  transverse  cracking  model  developed  earlier  in 
Section  II. 

In  addition,  laminates  of  the  fora  (Oj^O^/Oj]  were  tested  under  different 
temperature  environments. 

(c)  Transverse  Cracks  and  Delsminations  in  [+25/90n]s  -  Family  . 

This  series  of  tests  included  laminates  in  the  form  [+25/90nJa,  n  -  1/2,  1,2, 

3, 4, 6, 8.  The  reason  for  selecting  the  +25*-layer  as  the  major  load-carrying 

element  is  that  it  offers  the  optimal  lamina  constraint  which  Induces  free-edge 
delamination  under  uniaxial  tension.  Indeed,  as  will  be  shown  by  the  experi¬ 
mental  results,  the  laminates  of  n  -  1/2  and  1  suffered  free  edge  delamination 
without  any  transverse  cracking;  and,  the  delamination  mode  was  primarily  of 
mode-I.  This,  then  provides  an  ideal  case  study  for  the  delamination  modeling. 
As  for  those  laminates  with  n  >  1,  both  transverse  cracking  and  edge  delamina¬ 
tion  were  induced,  resulting  in  complicated  interacting  effects.  In  particular, 
for  cases  of  n  >  3,  other  modes  of  fracture  became  predominant.  Thus,  this 
series  of  tests  presented  a  multitude  of  Interesting  fracture  growth  processes. 
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In  addition  to  tha  [+25/90nlg  aariaa,  laminates  o£  vara  also 

tastad.  This  is  to  chack  whether  or  not  tha  ultimata  strength  of  tha  load 
carrying  layer  (tha  ^S^lmyar)  would  be  affected  by  tha  layer-thickness 
affect.  This  can  be  accomplished  by  comparing  it  to  the  [+25/90]  laminate, 
whose  load-carrying  +2 5- layer  has  only  half  the  thickness. 

(d)  Mixed  Mode  Fracture  in  Notched  [0g]  -  Family.  A  series  of  tests 
was  conducted  in  [0g],  ©  ■  0°,  15°,  30®,  45®,  60®,  75®  and  90®.  In  each  case 
the  specimen  had  double  side  notches  of  0. 15"  size.  Under  uniaxial  tension, 
fracture  failure  was  induced  along  the  fiber  direction.  The  fracture  is  of 
mixed-mode  in  general,  and  the  ratio  between  mode-I  and  mode-II  depends  on 
the  value  of  ©.  Results  of  this  test  has  provided  some  insights  about  the  mater' 
ial  resistance  under  various  degrees  of  mixed-mode  fracture. 

3.2  Experimental  Procedures 

Throughout  this  study,  all  laminates  tested  were  fabricated  from  Fiberlte 
T300/934  prepreg  tape.  They  were  autoclave  cured  at  350®F  (450°K)  according 
to  the  manufacturer's  specifications.  Nominal  fiber  volume. was  determined  in 
the  range  of  66  +  22.  The  12"xl2"  laminated  panels  were  cut  by  diamond  saw 
to  a  coupon  length  of  9"  and  width  of  1".  Fiberglass  tabs  of  length  1.5" 
were  bonded  to  each  end  of  the  coupon.  Coupons  were  stored  under  normal  room 
conditions  (70~80®F  and  60  +52  relative  humidity)  for  a  typical  duration  of 
3  to  6  months  before  test. 

All  tests  were  conducted  on  an  Instron  Universal  Tester  under  normal 


room  conditions  (except  otherwise  stated).  Load  was  applied  at  a  displace¬ 
ment-control  rate  of  0.01"  per  minute. 


The  experimental  procedure  for  documenting  daaege  sequence  In  the  various 


laminates  was  as  follows: 

(1)  The  ultimate  tensile  strength  (UTS)  for  each  laminate  in  the  series  was 
determined  by  a  minimum  of  four  ramp  tensile  tests.  These  specimens  were 
strain  gaged  to  determine  the  tensile  stiffness  properties,  including  the 
laminate's  Young's  modulus  and  Poisson  ratios. 

(2)  For  each  laminate  type,  a  minimum  of  four  specimens  were  step-loaded 
under  axial  tension.  The  initial  step  loading  Increment  was  approximately 
10Z  UTS  until  damage  in  the  form  of  transverse  cracks  and/or  de lamina¬ 
tion  was  detected.  Then  the  load  increment  was  subsequently  reduced  to 
approximately  5Z  UTS  until  the  ultimate  laminate  failure. 

(3)  After  each  step  loading,  the  specimen  was  removed  from  the  testing 
machine  and  was  inspected  for  damage  by  means  of  DIB  (dl-lodobutane) 
enhanced  x-radlography . 

(4)  In  addition  to  x-ray  examination,  the  free  edges  of  each  specimen  were 
also  optically  examined  using  a  microscope  to  observe  the  details  of  the 
damaged  area. 

(5)  In  many  cases,  tested  specimens  were  selectively  sectioned,  polished 
and  examined  under  SEM  (Scanning  Electronic  Microscope) .  The  SEM  exami¬ 
nation  provides  additional  internal  cracking  details  in  the  thickness 
sectional  view  of  the  laminates. 

Figure  3.1  Illustrates  a  typical  x-ray  sequence  taken  during  the  incre¬ 
mental  testing  process  for  a  [+25/903]g  specimen.  Here,  the  load-damage 
relationship  of  the  specimen  can  be  studied  in  detail  by  examining  the  series 
of  x-radiographs.  In  this  case,  damage  of  any  form  was  not  detected  up  to 
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Flgura  3.1  X-radiographa  Showing  Tranavarse  Cracking  and  Sdg«  Dalaainatlon 
Growth  Soluanoa  Undar  Stop-load  Fbr  A  C +25/90 0  Spaolaaa. 
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a  load  of  SOX  the  ultimate  tensile  strength  (UTS).  The  first  sign  of  damage 
was  a  transverse  crack  which  appeared  at  55X  UTS,  identified  by  the  thin  dark 
line  perpendicular  to  the  load  direction  in  the  x-ray  picture.  More  cracks 
were  seen  to  appear  in  the  next  several  load  increments. 

At  the  load  of  88X  UTS,  the  first  sign  of  edge  delamination  appeared  in 
the  form  of  a  thumbnail  shaped  crack  at  the  left  side  of  the  specimen.  The 
area  of  delamination  extended  stably  as  the  applied  load  was  Increased.  From 
the  x-ray  pictures,  it  is  seen  that  a  higher  density  of  transverse  cracks  exists 
in  the  delaminated  area.  And,  by  comparing  the  pictures  at  99Z  UTS  and  100Z 
UTS,  it  is  evident  that  the  final  failure  was  quite  sudden  and  it  did  not  seem 
to  involve  extensive  delamination  in  this  particular  sample  (the  reason  for 
this  growth  behavior  will  be  explained  later  in  this  section). 

Optical  edge-view  examination  of  the  delaminated  specimen  showed  that 
the  delamlnatlon  crack  in  this  case  followed  primarily  the  -25°/90°  interface. 

Fig.  3.2  shows  another  sequence  of  x-radiographs  (enlarged)  taken  for  a 
[+25/90] 8  specimen.  It  is  seen  here  that  edge  delamination  grew  much  more 
extensively  before  final  failure.  Furthermore,  in  this  case,  no  transverse 
cracks  were  detected  before  delamination.  And,  when  they  appeared,  they 
existed  only  in  the  delaminated  area.  The  delamlnatlon  crack  was  observed  to 
zig-zag  within  the  90a-layer  in  an  edge-view  optical  examination. 

From  che  load-sequence  x-ray  pictures,  either  the  transverse  crack 
density  or  the  percent  of  delamlnatlon  may  be  plotted  as  a  function  of  the 
applied  load.  For  example,  a  relationship  between  the  transverse  crack  density 
(number  of  transverse  cracks  per  unit  length,  N)  and  the  applied  load  a  can 
be  obtained  from  the  x-ray  sequence  shown  in  Figure  3.1.  This  relation  is 
shown  graphically  in  Figure  3.3.  The  load  required  to  initiate  transverse 


56 


Figure  3.2  X-radiographa  of  A  [±25/90l  Specimen  Taken 
and  (b)  100*  UTS.  * 
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Figura  3.3.  Tranavarae  Crack  Danaity  va.  Appliad  Tanaila  Straaa 
for  a  [±25/903)#  Spaciaan. 
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cracking  o^c  is  determined  by  the  o-intercept  of  the  curve  as  shown  in  the 
figure.  As  previously  discussed,  at  the  onset  of  free-edge  delamination  a 
higher  density  of  transverse  cracks  is  observed  in  the  delaminated  region  of 
the  laminate.  This  Indicates  that  there  is  a  certain  interaction  between 
the  two  cracking  events.  Tims,  the  measurement  of  the  N-o  relationship  is 
terminated  at  the  onset  of  edge  delamination. 

The  area  of  free-edge  delamination  can  also  be  related  to  the  applied 
load  from  these  experimental  results.  The  onset  load  for  free-edge  delamina¬ 
tion  may  then  be  obtained  in  a  manner  similar  to  the  method  that  determines  the 
onset  load  for  transverse-cracking  as  described  above. 

As  has  been  mentioned  earlier,  a  number  of  specimens  were  also  sectioned 
for  SEM  examination.  Typical  transverse  cracks  in  the  90°-layer  can  be. seen 
in  a  thickness  sectional  view,  such  as  shown  in  Fig.  1.3.  Similarly,  edge 
sectional  view  of  the  laminate  gives  the  internal  details  of  the  edge  crack, 
such  as  shown  in  Fig.  1.5. 

3.3  Experimental  Results. 

For  all  the  results  reported  in  this  section,  each  represents  the 
averaged  values  from  three  or  four  (sometimes  more  than  four)  replicate 
specimens.  Generally,  data  scatter  is  in  the  5%  (standard  deviation)  range 
for  the  compliance  properties,  10Z  range  for  the  strength  properties.  In  this 
section,  the  statistical  nature  of  the  data  scatter  will  not  be  discussed. 
However,  whenever  appropriate,  the  range  of  the  data  scatter  will  be  Indicated. 
Otherwise,  only  the  averaged  value  will  be  used. 

(a)  Basic  Plv  and  Laminate  Properties 

For  the  T300/934  graphite-epoxy  composite  material  used  in  this  study, 
a  routine  material  property  characterization  test  was  conducted.  Some 
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details  of  the  characterization  test  are  found  in  Ref.  [9].  We  shall  mention 
only  that  the  unidirectional  fiber  reinforced  system  was  considered  macro- 
scopically  a  homogeneous  medium  which  possesses  a  monocllnic  symmetry  in  its 
principal  material  coordinates  (L.T.z)  [27].  The  engineering  stiffness  con¬ 
stants  determined  using  8-ply  laminates  and  tested  under  nominal  room  con¬ 
ditions  are  as  follows: 


el- 

21.0  Msi 

tj  -1.70  Msi 

Ez 

-1.70  Msi 

VLT  ’ 

0.3 

VTZ  ‘  0,54 

VLZ 

-  0.3 

glt  “ 

0.94  Msi 

Gtz  -  0.5  Msi 

glz 

-  0.94  Msi 

aL  * 

0.2  ue/°F 

aT  -  16  ye/°F 

az 

-  16  ue/°F 

(3.1) 


Each  material  ply  when  cured  in  a  laminate  has  a  nominal  thickness, 

t  -  0.0052"  (0.0132  cm)  (3.2) 

The  strength  property  of  the  unidirectional  system  is  given  in  Table 
3.1.  Since  both  the  longitudinal  and  the  transverse  stress-strain  relations 
(under  uniaxial  tension  and  room  conditions)  are  essentially  linear,  except 
near  failure,  it  is  sometimes  more  convenient  to  use  the  ultimate  strain  eu 
to  characterize  the  strength  of  the  material.  Note,  however,  that  the 
physically  measured  quantity  is  always  the  far  field  applied  load;  the  term 
ultimate  strength  or  ultimate  strain  e^  («  <JU^X)  represents  only  an  average 
over  the  gross  dimension  of  the  test  specimen. 
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Table  3.1  Test  Data  of  Some  Basic  Laminates 


Laminate 

E^,  msi 

o  ,  ksi 
u 

V  * 

[o8J 

21.0(20  -  23.5) 

_  250(231  -  270) 

1. 1(1.0  - 

1.15) 

[90s] 

1.70  (1-65  -  1.77) 

6. 4(6.0  -  7.0) 

0.39(0.34  - 

0.44) 

l±«)2s 

11.1(11  -  11.15) 

75(70  -  78) 

0.67(0.66  - 

0.68) 

3.3(3.26  -  3.41) 

20(18.6  -  22) 

0.61(0.56  - 

0.67) 

3.15(3.04  -  3.27) 

18.1(17.5  -  19.5) 

0.57(0.55  - 

0.62) 

Results  from  testing  the  C^25] 2s*  and  l+45g]a  laminates  are  also 

summarized  in  Table  3.1. 

The  longitudinal  stress-strain  relation  of  the  [+25]^  specimen  maintained 
a  perfect  linearity  up  to  failure.  The  specimen  was  severed  in  the  mid- 
section  of  the  specimen,  involving  predominantly  fiber  breakage.  Fig.  3.4, 

It  should  be  noted  that  the  stacking  sequence  of  the  [+25] specimen  is 
actually  [25/-25/-25/25]  .  It  is  to  be  distinguished  from,  say  [+25- ]  which 
is  actually  1 252/ —2,52 1  a  •  The  failure  behavior  between  the  two  can  be 

drastically  different. 

The  reason  for  the  difference  in  the  failure  mechanisms  of  the  1+25] 
and  [+25.]  laminates  is  due  to  the  layer-thickness  effect.  This  same  effect 

—  i  8 

is  illustrated  more  vividly  by  the  strength  results  of  the  s  end  the 
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[+45.J a  laminates,  see  Table  3.1.  It  is  noted  that  the  average  strength  of 
the  former  Is  20  leal  while  the  average  strength  of  the  latter  Is  only  18  leal, 
a  reduction  of  10%  because  of  an  Increase  in  material  layer  thickness. 

A  closer  examination  of  the  failure  process  of  the  previous  two  lami¬ 
nates  indicates  that  failure  modes  in  both  cases  involved  primarily  matrix 
cracking,  not  fiber  breakage.  Fig.  3.5  shows  the  x-radiographs  of  a  s 

laminate  which  is  loaded  near  final  failure.  At  about  98%  UTS,  a  single 
surface  45°-crack  is  induced,  Fig.  3.5(a),  which  is  almost  immediately  followed 
by  multiple  -45* -clacks  in  the  inner  -45°-layer ,  Fig.  3.5(b).  It  may  thus 
be  concluded  that  the  formation  of  the  first  major  45°-crack  initiated  the 
process  of  failure  of  the  laminate.  Clearly,  the  formation  of  the  45*-crack 
must  be  a  fracture  phenomenon,  which  depends,  to  some  extent,  on  the  thickness 
of  the  layer  in  which  the  crack  is  formed. 

The  final  failure  mechanisms  in  these  basic  laminates,  however,  will 
be  Investigated  in  a  future  study.  They  are  outside  the  scope  of  the  present 
report. 

Finally,  experiments  on  the  residual-stress-free  temperature  yielded  a 
value  of  Tq  ranging  from  300°-3258F.  In  this  experiment,  an  unsymmetric 
[0^/90^]  laminate  was  used.  When  cured  at  350*F,  the  laminate  configuration 
was  flat.  After  cooling  down  to  room  temperature,  the  laaiinate  curved  into 
an  arc  as  shown  in  Fig.  3.6.  Thus,  upon  re-heating  the  laminate  at  a  higher 
temperature,  the  curvature  of  the  specimen  decreases.  By  optically  measuring 
the  laminate  cord  height,  h,  at  various  test  temperature,  it  is  then 
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Figure  3.5  Failure  Mode  of  a  [+A53] s  Specimen,  (a)  A  Major 
45#-Crack  Formed  at  98Z  UTS,  (b)  Multiple  -45*- 
C racks  Formed  at  99%  UTS. 


■ftt  if  v».  -  ^4  -fM *  V* 


Laminates 


possible  to  determine  experimentally  the  residual-stress-frae  temperature  TQ  . 
Fig.  3.5  shows  the  measured  cord  height  h  versus. .t  separatum;  the  test  was 
conducted  for  s  specimen  immedlstely  after  curing.  It  is  seen  that  a  loss 
of  about  10*F  is  resulted  due  to  sane  Initial  stress  relaxation.  Upon  storage 
in  room  conditions  for  3-6  months,  further  relaxation  due  to  time  and  moisture 
absorption  was  encountered.  This  caused  an  additional  decrease  in  the  stress- 
free  temperature.  For  simplicity,  and  without  causing  significant  discrepancy, 
the  stress-free  temperature  Tq  for  all  laminates  tested  in  this  report  is 
approximated  by  a  constant  value  of, 

Tq  -  300*F  (3.3) 

(b)  Crack  Growth  in  the  [0/90n/0]  -  Family  In  thi8  family,  the  laminates 
[0/90n/0],  n  ■  1.2. 3. 4  and  the  laminates  [O^/PO^/O^,] ,  n  "  w*r* 
under  room  conditions.  All,  except  the  [0/90/0]  laminates,  showed  transverse 
cracks  under  the  uniaxial  tensile  load.  In  the  case  of  the  [0/90/0]  specimens, 
no  transverse  cracking  was  observed  until  the  specimens  were  severed  at  final 
failure.  Table  3.2  summarizes  the  test  results,  including  the  stiffness  Ex» 
onset  stress  for  transverse  cracking  and  the  final  failure  stress  of  the 
various  laminates. 

Plots  of  transverse  crack  density  N  versus  the  applied  load  o  were  ob¬ 
tained  from  the  loed-sequence  x-radlogrsphs  of  all  specimens,  except  those  of 
[0/90/0]  lamination.  These  plots  are  shown  in  Fig.  3.7.  It  is  seen  that  for 
the  cases  of  n  ■  2  and  3,  N  Increases  monotonlcally  with  a  until  final  failure; 
there  is  no  evidence  of  the  existence  of  the  so-called  "characteristic"  spacing 
discussed  earlier  in  Sections  I  end  II.  However,  for  the  esse  of  n  •  4,  a 
slight  leveling  off  of  N  is  indicated.  But,  for  this  case,  the  theoretical 
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characteristic  spacing  Is  S~ 16  t,  and  the  corresponding  characteristic  crack 
density  would  be  N  -  1/S  •'■'12.  From  Fig.  3.7,  the  leveling  off  crack  density 
is  more  than  25. 

Table  3.2  Test  Results  of  the  [0/9 0^/0] -Family 


Laminate 

E^.,  msi 

T.C.  onset  stress, 
kai 

Final  failure  stress, 
ksl 

[0/90/0] 

13.9 

No  visible  cracks 

144(126-158) 

[0/902/0] 

11.3 

68(59.0-70.5) 

108(85-122) 

[0/903/0] 

8.9 

48(44-53) 

86(84-88) 

[0/904/0] 

7.9 

41(34-44) 

76(73.5-80.5) 

[02/902/02] 

13.8 

85C76-98) 

162(140-180) 

[02/904/02] 

11.1 

50(42-57) 

114(106-117) 

The  results  listed  in  Table  3.2  are  graphically  displayed  in  Fig.  3.8, 
and  Fig.  3.9.  In  particular,  compare  the  results  shown  in  Fig.  3.9,  and 
note  the  different  transverse  crack  onset  load  between  [0/90/0]  and  [Oj/SOj/O], 
and  between  [0/902/0]  and  [0^1 90^/0^] .  The  effect  of  the  90“ -layer  thickness 
on  the  transverse  crack  formation  is  clearly  illustrated.  The  final  failure 
of  the  laminates  is  basically  determined  by  the  load-carrying  0*-layer; 
and  the  thickness  of  the  90*-layer  does  not  seem  to  have  any  effect  on  final 
failure. 
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Comparison  of  Test  Results  between  [0/90^0]  and  [O2/' *va/02]  Laminates 


Finally,  a  number  of  the  [Oj/^O^/Oj]  laminates  were  tested  under  con¬ 
trolled  temperature  environments.  In  addition  to  tests  in  room  temperature 
(78*F) ,  tests  were  also  conducted  under  -80*F  and  156°F  conditions.  Test 
results  are  displayed  in  Fig.  3.10.  It  is  seen  that  the  onset  stress  for 
transverse  cracking  increases  with  the  test  temperature,  suggesting  that  a 
toughening  effect  against  90°-cracking  is  developed  by  the  increase  temperature. 
However,  these  results  are  only  exploratory  in  nature,  more  comprehensive  ex¬ 
periments  are  needed  to  draw  any  concrete  conclusion.  The  final  failure  of 
the  laminate  is  controlled  by  the  0°-  layer,  and  it  is  not  effected  by  the 
increase  of  the  test  temperature. 

(c)  Crack  Growth  in  the  [+25/90n]s-Family.  in  this  series,  laminates  of 
[+25/90  J  ,  n  -  1/2,  1,2, 3, 4, 6, 8  as  well  as  [±257/90_]  were  incrementally 
loaded  and  radiographically  examined.  Table  3.3  summarizes  the  results  ob¬ 
tained  in  these  tests.  For  purpose  of  comparison,  tests  results  for  [+25] 
and  [90gJ  laminates  are  also  listed  in  Table  3.3.  The  average  stresses  at 
onset  of  transverse  cracking,  edge  delamination  and  final  failure  are  tabu¬ 
lated,  with  the  range  of  scatter  given  in  parenthesis. 

In  all  cases,  the  laminate  stress-strain  response  was  found  to  be 
essentially  linear  until  failure.  It  is  thus  convenient  to  use  the  laminate 
tensile  strain  as  a  measure  of  the  applied  load.  A  graphical  representa¬ 
tion  of  the  data  given  in  Table  3.3  is  displayed  in  Figi  3.11.  Here,  the 
applied  laminate  strain  at  the  onset  of  transverse  cracking,  edge  delamina¬ 
tion  and  final  failure  is  plotted  against  the  number  of  90° -plies,  denoted  by  n. 
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(90d)  Unidirectional  Laminate 


Figure  3.11  Observed  Tensile  Strain  at  the  Onset  of  Transverse  Cracking  and  Free' 
Edge  Delanlnation,  and  at  the  Final  Failure  vs.  the  Nuaber  of 
90*-plies  for  the  [+25/90n]8  Laalnates. 


Both  the  averaged  value  and  the  range  of  the  scatter  are  shown  for  each  case. 
The  strain  at  failure  Is  seen  to  be  approximately  constant  for  the  cases  of 
n  *  1/2,  1,2,3.  This  value  is  roughly  equal  to  the  tensile  strain  at  failure 
of  the  [+25]  2g  laminate.  For  larger  n,  the  strain  at  failure  drops  off 
significantly  (discussions  on  this  point  are  presented  in  Sections  IV  and  V). 

The  onset  strain  for  edge  delamination  also  appears  constant  for  the 
n  -  1/2,  1,2,3  laminates  and  it  drops  off  significantly  for  n  >  3.  The 
free-edge  delamination  generally  initiates  at  about  90Z  of  the  strain-at- 
fallure.  In  fact,  the  onset  of  edge  delamination  and  the  final  failure  were 
almost  simultaneous  in  the  [+25/90.]  laminates. 

—  4S 

Transverse  cracks  were  not  observed  in  the  n  *»  1/2  and  n  ■  1  laminates 
prior  to  the  onset  of  edge  delamlnatlon;  and  the  transverse  cracks  were 
largely  confined  within  the  delaminated  area  when  they  occurred.  In  the 
cases  of  n  2 ,  transverse  cracks  appeared  before  edge  delamlnatlon.  The 
onset  strain  for  transverse  cracking  decreases  in  laminates  of  n  *  2,3,4; 
and  then  Increases  slightly  for  n  »  6  and  n  ■  8.  Generally,  for  the  cases 
of  n  >  2,  the  onset  strain  for  transverse  cracking  is  less  than  the  ultimate 
tensile  strain  of  the  [9Qa]  laminate.  It  will  be  shown  in  Section  IV  that 
this  is  caused  by  the  residual  thermal  curing  stresses  in  the  laminate. 
Clearly,  the  onset  strain  for  transverse  cracking  shows  a  strong  dependence 
on  n. 

It  is  also  noted  that  the  onset  strains  for  transverse  cracking, 
edge  delamlnatlon,  and  the  final  failure  for  the  [+252/-252/902]8  laminate 
are  below  the  values  observed  for  the  [+25/90]  laminate  even  though  these 

9 

two  laminates  are  essentially  of  the  same  construction  except  for  the  layer 
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thickness.  Again,  the  layer  thickness  in  this  laminate  plays  an  Important 
role  in  both  the  cracking  processes  and  the  final  failure  mechanisms. 

In  fact,  the  strains  at  the  onset  of  transverse  cracking  for  the  [+252 /-252 /902 ] g 

and  the  [+25/-25/902] ^  laminates  are  about  the  same.  Interestingly,  the 
thickness  of  the  90* -layer  is  the  same  in  both  cases. 

The  transverse  crack  density  as  a  function  of  the  applied  load  is 
summarized  in  Figure  3.12  for  the  cases  of  n  ■  2, 3, 4, 6.  The  range  of  the 
data  is  shown  by  the  shaded  band  for  each  case.  The  transverse  crack  density 
N  is  seen  to  rise  sharply  following  the  onset  of  transverse  cracking.  The 
process  of  transverse  cracking  is  enhanced,  however,  by  free-edge  delamination. 
Thus,  an  experimental  limit  for  the  transverse  crack  density  could  not  be 
defined  in  these  laminates  either. 

'  With  this  overview,  it  is  seen  that  the  development  and  transition  cf 
damage  modes  is  dependent  upon  thickness  of  the  90*-layer.  It  is  now  con¬ 
venient  to  separate  the  laminates  of  the  [+25/90  ]  series  into  three  groups 

**  n  8 

as:  (1)  thin  90*-layer  (n  ■  1/2,  1);  (2)  thick  90°-layer  (n  -  2,3);  and  (3) 
very  thick  90*-layer  (n  •  4,6  8).  The  damage  development  results  for  these 
three  groups  will  be  discussed  in  more  detail  next. 

Thin  90°-layer  Laminates  (n  ■  1/2,  1) 

The  thin  90*-layer  laminates  did  not  suffer  transverse  cracks  prior 
to  delamination.  Their  edge  delamination  occurred  primarily  Inside  the 
90*-layer  of  the  laminate;  and  it  propagated  stably.  Final  failure  occurred 
at  the  laminate  strain  approximately  equal  to  the  failure  strain  of  the 
[±2532a  laminates. 
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Figure  3.12  Density  of  Transverse  Cracks  vs.  Applied  Tensile 
Stress  In  [+25/90n]s  Laminates  for  n  ■  2, 3, A  and  6 


The  damage  development  in  the  [+25/90 1/2  ]  laminates  (n  *  1/2) 
started  with  a  Tree-edge  delamination  along  the  length  of  one  side  of  the 
test  coupon  and  then  grew  larger  with  delaal nation  on  both  sides.  X-radlo- 
graphs  of  three  specimens  loaded  to  an  axial  strain  of  %x  -  0.63Z  are  shown 
in  Figure  3.13.  These  x-rays  show  what  appears  to  be  transverse  cracks  in 
the  delaminated  region.  As  the  ultimate  load  was  approached,  de laminations 
on  both  sides  were  observed  to  grow  stably  until  they  coalesced  in  the  middle 
of  the  test  coupon  and  resulted  in  final  failure.  This  is  shown  by  the  failed 
specimen  in  Figure  3.13.  An  oblique  section  cut  from  thin  failed  specimen 
was  SEM  examined.  A  magnified  section  view  of  the  SEM  picture  shows  that 
the  delamlnatlon  crack  seen  in  the  x-radiograph  actually  occurred  mainly  Inside 
the  90°-ply,  see  Fig.  3.14.  No  transverse  cracks  were  found  to  form  In  this 
laminate  even  after  failure.  The  delamlnatlon  is  seen  to  Initiate  In  the 
mid-plane  at  the  free-edget  but  grew  Inward  along  the  25/90  interface.  Several 
oblique  cracks  are  seen  in  the  90°-ply  which  may  be  responsible  for  the 
delamlnatlon  shift.  The  orientation  of  the  oblique  cracks  strongly  suggest 
that  a  locally  complex  state  of  stress  exists  In  the  vicinity  of  the  25/90 
Interface  delamlnatlon  crack  tip. 

In  the  case  of  the  [+25/90]  (n  ■  1)  laminate,  mid-plane  delamlnatlon  was 
also  observed  as  Che  first  form  of  damage.  However,  transverse  cracks  occurred 

In  the  delaminated  region,  see  Fig.  3*15 1  and  these  extended  across  the  entire 
width  of  the  specimen.  The  strains  at  the  onset  of  delamlnatlon  and  at  the 
final  failure  are  approximately  equal  to  those  for  the  n  •  1/2  case.  It  appears 
that  the  transverse  cracks  In  the  [+25/90]  laminate  did  not  have  any  effect 
on  the  mechanisms  that  caused  edge  delamlnatlon.  And  the  final  failure  Is 
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(b) 


Figure  3.14  Photomicrograph  of  tha  Delaminated  Region  in  the  [±25/90i^2]e 
Specimen  for  the  Oblique  Section  Shown  in  Fig.  3.13. 
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controlled  by  the  load-carrying  +25° -layer .  Tha  sub-laminate  cracking 

event  did  not  eeea  to  af£act  the  load-carrying  capacity  of  the  25°-layer . 

Thick  90*-laver  (n  ■  2.3) 

The  thickneae  of  the  90*-layer  In  theee  laninates  la  typically  found  in 
■oat  practical  laminates.  For  both  caaaa  of  n  •  2  end  n  »  3,  traneverae  crack¬ 
ing  was  the  flrat  damage  node  under  the  load.  A  high  denaity  of  trane¬ 
verae  cracks  vaa  reached  prior  to  the  onset  of  delaalnatlon.  Delamination 
occurred  mainly  in  the  90* -layer  for  n  -  2,  but  mainly  along  the  25/90  inter¬ 
face  for  n  -  3.  Delamination  crack  growth  was  generally  stable,  although 
the  growth  was  quite  limited  in  the  case  of  n  -  3,  see,  e.g..  Fig.  3.1.  The 
ultimate  atrength  of  the  n  ■  2  laminates  la  also  approximately  equal  to  the 
strength  of  the  [+25] 2>  laninatee.  But,  &  -reduction  In  the  final  failure 
strain  la  Indicated  for  the  case  of  n  »  3,  see  Fig.  3.11.  A  degradation 
In  strength  has  also  been  found  in  the  laminates  to  be  discussed  next. 

Very  Thick  90* -layer  Laninates  (n  ■  4.6.8) 

The  first  form  of  damage  for  the  very  thick  90* -layer  laminates  is  also 
transverse  cracks,  which  occur  at  a  much  lower  stress  level.  Stress  con¬ 
centrations  at  the  transverse  crack  tip  and  at  the  frae-edge  of  the  specimen 
cause  a  delaalnatlon  along  the  25/90  interface.  The  triangular  shape  of  the 
delaalnatlon  plans  observed  in  the  x-radlo graphs  indicates  that  there  is  an 
interaction  between  these  two  modes  of  damage.  The  growth  of  this  type  of 
delaalnatlon  appeared  very  unstable;  and  the  final  failure  occurred  at  a 
laminate  strength  significantly  lass  than  the  strength  of  the  [+25] laminate 
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The  delamination  growth  process  of  a  [+25/90^]  laminate  Is  illustrated 
by  the  x-radiographs  shown  in  Figure  3. 16.  The  x-rays  show  that  no  typical 
edge  delaalnation  occurred  up  to  a  load  of  96. 8Z  UTS.  At  the  load  of  99. 91  UTS, 
a  delaalnation  originated  at  the  transverse  crack  tips  is  developed,  which 
extends  across  the  full  width  of  the  coupon.  A  small  increase  in  load  was  then 
sufficient  to  cause  the  final  failure.  An  important  feature  of  the  damage  in 
this  (n  •  4)  laminate  is  the  triangular  shape  of  the  delaalnation  near  the 
free  edge.  The  apex  of  the  triangle  is  clearly  aligned  with  the  transverse 
crack  shown  by  the  arrow  in  the  figure.  In  the  previously  discussed  cases 
of  n  <  4,  the  delaalnation  crack  front  was  always  of  a  thumbnail  shape  when 
it  first  formed.  This  difference  in  the  shape  of  delaalnation  indicates 
clearly  a  change  in  the  mechanisms  that  control  the  initiation  and  growth  of 
the  delaalnation.  The  damage  development  process  for  the  case  of  n  •  6  is 
similar,  with  more  pronounced  transverse  crack-induced  delaalnation  which  was 
observed  in  the  experiment  of  this  laminate. 

Figure  3.17  displays  the  daaage  In  two  [+25/90g]s  laminates  after 
tensile  loading  to  97X  UTS.  An  example  of  transverse  crack  tip  induced 
delamination  is  shown  in  picture  (a),  where  a  dark  sons  (delaalnation) 
associated  with  a  transverse  crack  is  Indicated  by  an  arrow.  The  x-radio- 
graph  in  (b)  shows  delaalnations  occurred  in  both  the  upper  and  lower  25/90 
interfaces.  The  apex  of  the  triangular  shaped  dark  areas  marked  by  the 
arrows  are  clearly  aligned  with  a  transverse  crack  tip  delamination.  It  is 
concluded  that  the  stress  concentrations  at  a  transverse  crack  tip  interact 
with  the  free-edge  stresses  to  cause  this  delaalnation. 
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Figure  3.16  X --radiographs  Showing  Edge  Delaalaation 

Growth  in  a  [+25/904]g  Specimen.  \ 
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Figure  3.17  X-radiographs  of  Two  [+25/90gls  Specimens  after  Tenaile  Loading  to  97X  UTS 


the  foregoing  discussions  illustrated  the  significant  effect  of  the 

90* -layer  on  the  damage  modes  and  their  growth  nechanlsas  in  the  [+25/90  ] 

— ■  n  8 

laminates.  Fig.  3.18  provides  a  schematic  summary  of  the  fracture  sequences 
in  this  [+25/90  1  series.  Column  <a)  shows  the  development  of  transverse 
cracks  just  before  the  onset  of  edge  delamination,  as  seen  in  the  plane  view 
x-radlographa .  Column  (b)  shows  the  plane  and  edge  views  (the  broken  lines 
are  the  25/90  interface)  of  the  test  coupon  isanediately  after  edge  delamina¬ 
tion.  Column  (c)  shows  the  plane  views  and  edge  views  of  the  specimen  before 
their  final  failure  (transverse  cracks  are  shown  in  the  edge  views  only  when 
they  are  associated  with  the  initiation  of  a  25/90  delamination  at  the  trans¬ 
verse  crack  tip). 

The  fracture  sequences  shown  in  Fig.  3.18  illustrate  that  the  chosen 
grouping  by  90*-ldyer  thickness  coincides  with  transitions  in  the  several 
different  damage  modes  in  the  [+25/90Q]g  series  found  in  the  experiment. 

Recall  for  the  [£  25/90n]g  laminates,  that  the  observed  fracture  modes 
consisted  of  mainly  fiber  breakage.  The  failure  section  of  the  specimen 
was  essentially  transverse  to  the  load  direction,  see  Fig.  3.4.  Similar 
fracture  modes  were  observed  for  the  +25 '-layer  in  the  final  failure  of  all 
the  laminates  in  the  [+25/90^] ^  series.  The  morphology  of  the  fracture  of 
the  +2 5 '-layers  appears  to  be  the  seme  for  all  of  the  laminates  investigated. 
However,  the  laminate  strains  Sx  at  the  final  failure  show  a  strong  dependence 
on  the  thickness  of  the  90* -layer,  especially  when  it  is  thick.  Further, 
once  edge  delamination  occurs,  final  failure  is  imminent.  Apparently,  the 
growth  of  delamlnatlon  has  some  ill-effects  on  the  final  failure;  the  exact 
mechanisms  will  not  be  investigated  in  this  report,  however. 
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Figure  3.18  Schematic  of  the  Fracture  Sequences  in  the  [+25/90^] 
Laminates,  (a)  Just  Prior  to  Edge  Delamination, 

0>)  Subsequent  to  Edge  Delamination,  and  (c)  Just 
Prior  to  the  Final  Failure. 


(d)  Mixed-Mode  Fracture  in  Notched  [0g] -Family.  This  series  of  tests  con¬ 

sisted  of  a  class  of  off-axis,  [0g]  laminates.  The  test  specimens  were  cut 
with  0.15"  double-side  notches  using  a  nominal  8-mil  diamond  blade.  The  off- 
axis  angle  0  varied  from  0°  to  90"  at  15°  intervals.  A  general  configuration 
of  the  test  specimen  is  illustrated  in  Fig.  3.19 . 

Under  uniaxial  tension,  fracture  failure  was  induced  in  each  specimen 
with  the  crack  initiating  from  the  notch  tip  and  propagating  along  the  fiber 
direction.  Thus,  a  class  of  matrix-dominated,  mixed-mode  fracture  was  ex¬ 
perimentally  determined.  Fig.  3.19  shows  a  plot  of  the  far-fleld  laminate 
stress  Oq  at  onset  of  fracture  versus  the  fiber  angle  0. 

A  finite  element  simulation  of  the  crack  propagating  from  the  notch 
tip  along  the  fiber  direction  was  then  conducted  (the  numerical  details  are 
presented  in  Section  IV).  The  calculated  energy  release  rates  G^.  and  GII 
for  each  value  of  0  are  displayed  in  Fig.  3.20.  It  is  seen  that  for  0  _>  60 
the  fracture  is  essentially  of  mode  I  only.  For  smaller  0,  mode-II  contri¬ 
bution  increases  rapidly.  Fig.  3.21  shows  the  ratio  of  G^/G^  as  a  function 
of  0. 

Correlating  the  calculated  energy  release  rates  with  the  experimental 
results,  the  materials  fracture  resistance  property  may  be  evaluated.  For 
example,  Fig.  3.22  shows  the  experimentally  determined  mixed-mode  energy 
release  rate  G^  of  the  material  as  a  function  of  the  G^/Gj  ratio.  It 

is  seen  that  for  0  >.60°,  xDc^^Ic  bec8u8e  the  fracture  was  essentially 
of  mode  I,  the  average  value  for  GIc  is 

GIc^0.9  in-lb/in2  (3.4) 
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Figure  3.19  Mixed-Mode  Fracture  Data  of  the  Notched  [0g]  Laminates 
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Figure  3.22  Material  Fracture  Resistance  G(I>Ii)c  v8r8U8  Gn^Gi  Ratio 


Not*  that  the  value  of  G^c  determined  In  thia  experiment  rapraaanta 
the  fraetura  resistance  of  the  unidirectional  composite  against  the  type 
of  cracking  as  illustrated  in  Fig.  1.8(a);  it  is  essentially  a  case  of  0*/0* 
Interface  delawlnatlon. 

Similar  values  for  GIc  (0.8-0.92  in-lb/in2)  were  obtained  by  Cullen  [11] 
and  Williams  [12]  who  used  the  AS-3502  composite  system  tested  undet 
a  different  set  of  specimen  configurations. 

From  Fig.  3.22,  it  is  also  seen  that  the  value  of  G^  gl)c  increases 
with  the  ratio  of  G^/G^,  suggesting  that  the  energy  required  to  fracture 
a  mixed-mode  crack  is  not  constant;  it  may  be  dependent  on  the  amount  of  the 
shearing  action.  In  addition,  the  plot  in  Fig.  3.22  indicates  the  existence 
of  GJIc,  being 


GIIc~2.80  in-lb/in 


2 


(3.5) 


Dependence  of  the  fracture  energy  on  the  G^/Gj  ratio  has  been  encountered 
recently  by  Vanderkley  [13],  who  also  tested  the  unidirectional  AS-3502 
composite  system.  And,  the  existence  of  G^c  (at  the  macroscopic  level) 
has  been  discussed  by  Wilkins  [28].  For  the  AS-3501-06  graphite-epoxy 
system  used,  Wilkins  found  the  value  of  G^  ^^(OVO*  delamination)  ranged 
from  2.1  to  2.98  in-lb/in2  with  Gjj/ G^j.  equaled  3.4.  This  range  is  within 
the  same  experimental  range  of  the  present  study. 

Apparently,  the  problem  of  mixed -mode  crack  propagation  is  much  more 


complicated  than  the  limited  study  has  shown  here.  It  requires  concentrated 
efforts  on  this  particular  subject,  both  experimentally  and  analytically. 


In  chc  e««c  of  the  double-aide  notched  [90g]  laminate,  the  calculated 
energy  releaae  rate  Gj  can  be  expressed  in  terms  of  the  notch  else  a, 

GjU)  -  4.75xl06-a.Jx*  (3.6) 

Now,  using  the  failure  strain  of  the  unnotched  [90g]  epeclaen  ix“  0.39X 
froa  Table  3.1,  and  the  critical  energy  release  rate  Gjc=  0.9  lb/l11  f*°®  (3*^)» 
we  obtained  for  the  macroscopic  material  flaw  size 

aQ~  0.0125".  (3.7) 

This  Is  about  2.5  times  the  thickness  of  a  single  material  ply. 


IV.  FINITE  ELEMENT  ANALYSIS  AMD  EXPERIMENTAL  COgRELATggS 

4.1  Finite  Element  Fora^i  *«»<'»« 

Shown  In  Fig.  4.1  In  the  cross-section  of  a  symmetric  laminate.  4mm— 
tha  laminate  la  long  la  tha  axial  x-dir  action.  Than*  under  tha  axial 
tanalla  load  F^,  tha  dlaplaeaaant  flald  of  tha  laminate  can  be  expreaaed  In 
the  general  form  [3]. 


-  •fft  +  U(y,«) 

v(*.F.a)  ■  V(y,t)  (4.1) 

v(x,y,s)  -  W(y,a) 


where  e^  la  the  uniform  axial  laminate  a train  Induced  by  F^,  and  V  and  W 
are  the  y  and  z  dlaplacaaanta  respectively,  which  are  Independant  of  x. 

The  function  U  represents  the  warping  of  the  laminates  cross-section.  Owing 
to  synmatry,  the  solution  field  can  be  reduced  to  the  first  quadrant  of 
the  laminate  cross  section  as  shown  by  the  shaded  region  in  Fig.  4.1.  The 
solution  region  is  then  represented  by  a  network  of  nodes  and  triangular 
elements  as  shown  in  Fig.  4.2  with  the  following  symmetric  boundary  con¬ 
ditions. 


u(0,z)  -  v(0,z)  -  0 
w(y.O)  -  0 


(4.2) 


Tha  formulation  of  the  (constant  strain)  triangular  finite  element 
based  on  tha  displacement  field  of  (4.1)  begins  with  the  nodal  force-dis¬ 
placement  relationship  for  a  typical  element. 


B  •  B  v 
*1  *1 


Kij 


l,j  -  1,9 


(4.3) 
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where  K  la  tha  element  stiffness  matrix  and  F®  rapraaaata  tha  el  am  ant 
residual  nodal  force  vector.  Tha  6  vector  denotea  the  u,  v,  v  functlona 
at  the  three  nodea  of  the  element  and  R  repreaenta  the  "reactive**  nodal 
forcea  (see  detalle  in  Baf.  [6J). 

Assume  the  material  of  the  element  la  described  by  the  generalized 
Hooke* a  Lav  as 

'i  •  Vi  -  VT>  l-J  ■  1*‘  <*•*> 

Then,  the  element  stiffness  matrix  K  la  given  by  [6] , 


(BiA»Vv 


l. 1  -  1.9 

m, n  «  1,5 


and  the  residual  nodal  forces  R®  are  given  by 


(4.5) 


R° 


<Bi.D«.*n*  +  ‘in0. 


')T 


m,n 

1 


1.5 

1,9 


(4.6) 


where  B  la  a  5  x  9  operator  matrix  depending  on  the  coordinates  of  the 
three  nodes.  In  (4.5)  and  (4.6),  V  is  the  volume  (area)  of  the  element  and 


Dmn  "  C(mfl)(n+1) 

en*  -  -  a^jdT  m.n  «  1,5  (4.7) 

°n*  ’  <Wl)‘*0  ’  “l411 

If  the  axial  strain  Sq  is  known,  then  the  "residual"  nodal  forces  R° 
can  be  calculated.  Thus,  for  a  total  of  N  elements  in  the  network,  N  seta 
of  nine  equations  are  obtained.  However,  these  9N  equations  can  be  con¬ 
solidated  into  a  set  of  3M  linear  equations  In  the  fora 


fi  -  Vj  311 


(4.8) 
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H  being  the  nunber  of  nodes  in  the  network.  In  (4.8),  f  denotes  the  system 
resultant  nodel  fores  vector  which  consists  of  both  the  residual  force  R 9 
and  the  reactive  force  R;  A  denotes  the  global  stiffness  matrix;  and  d 
represents  the  displacement  vector  of  the  M  nodes  in  the  fields.  Equation 
(4.8)  can  be  solved  either  given  a  force  boundary  condition,  a  displacement 
boundary  condition,  or  nixed  boundary  conditions  for  the  nodal  network. 

The  local  stresses  and  strains  for  each  of  the  N  eleaents  can  then  be 
calculated  from  the  nodal  dlsplaceaent  solution  as 


*1  *  *0 

Vl  -  B-ld1  .  -  1,5;  1  -  1,» 

“i  '  CU<*j  -  V”  i  -  1.6 

Vl  "  »„<•.*  +  Vl>  +  *«*  1.5 


(4.9) 

(.10) 


If  the  applied  axial  force  is  specified,  then  the  value  of  e^  oust 
be  determined  to  satisfy  the  equilibriua  condition  that  the  sun  of  the 
forces  in  the  x-direction  be  equal  to  F^.  For  the  case  of  a  thermal  load 
condition,  the  applied  axial  force  F^  is  equal  to  rero.  Then,  eQ  must  be 
determined ;  eQ  represents  the  t hemal  axial  strain  in  the  plate  ^caused 

by  the  constant  thermal  load  AT. 

The  present  formulation  can  be  applied  to  solve  the  free-edge  de lami¬ 
nation  sad  transverse  cracking  processes  by  naans  of  the  "crack-closure" 
technique  described  in  Section  II. 


98 


The  delaminatioh  process  i«  modeled  by  taking  the  axial  x-dlreetfbu  as 
the  load  direction  of  the  laminate,  and  the  applied  laminate  load  la  C,  -  The 
free-* edge  of  the  laminate  is  at  y  •  L  as  shown  In  Fig.  4.2.  Then  the  delanlns- 
tlon  process  may  be  modeled  by  a  crack  starting  at  the  free-edge  and  extend¬ 
ing  towards  the  z-axls  along  a  chosen  plane  defined  by  z  •  constant. 

The  transverse  cracking  process  Is  modeled  by  taking  the  y-direction  as 
the  load  direction  of  the  laminate.  The  laminate  is  loaded  by  prescribing  a 
uniform  y-displacement  boundary  condition  along  y  »  L.  The  applied  load  in 
the  x-dir action,  F^,  is  specified  as  zero  to  maintain  a  uniaxial  tensile  load 
condition  in  the  y-dir action.  Then,  the  transverse  cracking  process  may 
be  modeled  by  a  crack  starting  at  z  *  0  and  extending  in  the  positive  z 
direction. 

Some  simplifications  in  the  calculation  of  G  may  be  made  when  a  crack 
is  extended  along  a  line  of  symetry.  The  details  concerning  these  sim¬ 
plifications,  the  loading  conditions,  and  the  boundary  conditions  will  be 
discussed  next  in  the  actual  calculation  of  the  transverse  cracking  and 
delamination  processes. 

4.2  Calculation,  Of  The  Energy.  Release  Rate 

In  the  present  analysis,  both  mechanical  and  thermal  loadings  are 
considered.  The  mechanical  load  is  the  uniaxial  tension  and  the  thermal 
loading  is  the  uniform  temperature  drop  AT,  Within  the  context  of  linear 
elasticity,  the  mechanical  and  thermal  loading  cases  may  be  solved  indepen¬ 
dently  under  unit  load  conditions  and  then  superimposed  to  obtain  the 
solution  for  a  combined  load  condition.  It  is  convenient  to  choose  the 
mechanical  load  as  the  far  field  laminate  strain  eQ  equal  to  one  microstrain 


(1  x  10  6).  This  represents  s  "fixed-grip"  load  condition.  The  theraal 
load  AT  is  chosen  as  minus  one  degree  Fahrenheit  for  convenience.  Let  f 

a 

end  be  the  nodal  forces  and  displacements  respectively,  due  to  the  unit 
mechanical  strain  load  as  determined  by  the  finite  element  solution  of 
Equation  (4.8).  Similarly,  let  7^  and  3^  be  the  nodal  forces  and  displace¬ 
ments  respectively,  due  to  the  unit  thermal  load  of  AT  -  -1°F.  Then,  under 
a  combined  loading  of  eQ  and  AT,  the  nodal  forces  and  displacements  can  be 
expressed  as 


f  -  £  e  +  f-AT 
m  l 

i  -  d  e  +  d_AT 
in  x 


(4.11) 


To  calculate  the  strain  energy  release  rate  It  Is  necessary  to  intro¬ 
duce  a  crack  in  the  finite  element  network  as  described  in  Section  II. 
Referring  to  Fig.  2.4,  let  D  be  the  relative  displacements  of  nodes  f  and  g 
when  the  crack  is  extended  to  a  +  Aa  defined  as 

f)  »  u-  -  u  (4.12) 

ml 

D  may  be  represented  In  terms  of  a  combined  mechanical  and  thermal  load 
of  e^  and  AT  as 

D  ■  !>  e  +  5-AT  (4.13) 

in  l 

mm  mm 

where  and  D^,  are  the  relative  displacements  associated  with  the  unit 
mechanical  and  unit  thermal  loadings  respectively.  Let  7  be  the  nodal 
forces  required  to  close  nodes  f  and  g  together.  F  may  be  expressed  in 
terms  of  a  combined  load  of  eg  and  AT  as 

f  -  fBe  +  f^T  (4.14) 
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where  F  and  F_  are  the  nodal  closing  £orces  associated  with  the  unit 

ID  1 

mechanical  and  unit  thermal  loadings  respectively.  The  energy  release 
rate  per  unit  thickness  may  then  be  expressed  as 

G  ■  fer  (v + VT)  •  <v +  <4-15) 

where  Aa  is  the  incremental  crack  extension  as  shown  in  Fig.  2.4. 

Equation  (4.15)  may  be  rewritten  in  the  form 

G-  [Ce(eQ)2  +  CT(AT)2  +  CeX(e)(AT)]t  (4.16) 

A 

where 

C  -  (f  •  5  )/2Aa 
e  mm 

CT  -  (Ft  •  DT)/2Aa  (4.17) 

CeT  “  (fm  *  5T  +  fT  ‘  V/2Aa 

The  parameter  t,  appearing  in  Eq.  (4.16)  is  the  scale  factor  of  length 
between  the  actual  dimension  and  the  finite  element  dimension  of  the  specimen. 
Note  that,  for  a  given  geometry  and  material  properties,  equations  (4.17) 
are  functions  of  the  non-dimensional  crack  length  a/t  and  are  independent 
of  the  applied  loading.  Thus,  the  functions  Ce,  C^,  and  C£X  are  character¬ 
istic  for  a  particular  crack  growth  and  are  hereafter  denoted  as  the  "energy 
release  rate  coefficient  functions". 

The  mechanical  load  coefficient  function  Cg  may  be  expressed  in  terms 
of  the  x,  y,  z  components  of  F  and  D  as 

ID  ul 

*”  The  expression  for  CeT  is  equivalent  to  Eq.  (2,31) 


D  )/2Aa 

my 

D  )/2Aa 


(4.18) 


noting  that 


C.I  *  <F.y 

°.II  ■  (P«x 


C.III  *  <p«*  ’  D«.)/aA* 


c  -  c  T  +  c 
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+  C 
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(4.19) 


The  subscripts  I,  II,  and  III  are  Identified  with  the  mode  1,  mode  II 
and  mode  III  crack  extension  modes  respectively.  Similar  expressions  may  be 
obtained  for  the  thermal  and  mixed  load  coefficient  functions.  Thus,  the 
mode  I,  mode  II,  and  mode  III  components  of  the  energy  release  rate  may  be 
expressed  in  terms  of  the  coefficient  function  as 

Gx  -  t[CeI(e)2  +  Cti(AT)2  +  Cen(e)(AT)] 

GII  "  tICeII(e)2  +  CTII(AT)2  +  GeTII(e)CAT)1  (4*20) 

°XII  -  t[c.m(e)2  +  ctiii(at)2  + 


4.3  Numerical  Solution  Convergence  Study 

The  numerical  solution  of  the  transverse  cracking  and  delamination 
problems  requires  an  optimal  finite  element  discretization  which  yields 
sufficiently  accurate  results  with  minimum  computational  requirements. 
Lacking  an  exact  solution  for  comparison,  the  acceptability  of  the  numerical 
approximation  is  determined  here  by  comparisons  between  the  results  obtained 
from  various  node  and  mesh  distributions. 


For  example,  consider  the  Interlaminar  shear  stress  distribution  in 
the  vicinity  of  a  transverse  crack  of  the  [+25/90]  laminate,  see  Fig.  2.10 
The  calculated  interlaminar  shear  stress  t  along  the  25/90  interface  is 

yz 

plotted  against  the  logarithm  of  the  distance  from  the  transverse  crack  tip 
in  Fig.  4.3.  It  is  seen  that  by  using  different  crack  tip  nodal  spacing*  the 
convergence  of  the  shear  stress  near  the  crack  tip  can  be  evaluated.  For  a 
spacing  of  da  ■  0.025t  (t  being  a  single  ply  thickness)*  the  semi-log  plot 
of  the  shear  stress  x  is  almost  a  straight  line;  the  slope  of  the  line  is 
nearly  -0.5.  This  agrees  with  the  commonly  accepted  notion  that  the  stress 
singularity  at  a  crack  tip  in  an  isotropic  medium  is  associated  with 

Similarly*  the  calculated  energy  release  rates  for  a  transverse  crack 
growing  in  the  [+25/90  ]  laminates  are  plotted  against  the  increment  crack 
length  Aa/nt  in  Fig.  4.4.  Here,  it  is  seen  that  for  Aa  -  0.05t*  the  calcu¬ 
lated  Gj  reaches  almost  99%  the  converged  value  (using  Aa  ■  0.025t). 

Apparently,  the  numerical  convergence  of  G  is  faster  than  the  crack-tip  stress* 
Fig.  4.5.  This  is  due  to  the  fact  that  G  is  calculated  from  the  nodal  forces 
and  displacements,  while  the  stress  is  calculated  from  the  first  derivatives 
of  the  nodal  displacements. 

Use  of  a  nodal  spacing  Aa  ■  0.050t  is  sufficient  for  purpose  of  ex¬ 
perimental  correlation  of  the  fracture  processes  discussed  in  this  study.  It 
should  be  noted  that,  the  single  ply  thickness  of  the  material  is  nominally 
t  -  0.0052",  Eq.  (3.2).  Then,  Aa(-  0.05t)  is  in  the  order  of  250  x  10-6  inch, 
which  is  about  2-3  times  the  graphite  fiber  diameter.  Tims*  the  choice  of  Aa 
is  probably  the  limit  for  which  the  assumption  of  "ply- elasticity"  remains 
valid. 
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Figure  4.4  Energy  Release  Rate  vs.  Crack  Length  for  Transverse  Crack  Mesh  Analysis 


4.4  Correlation  with  the  [ 0/90^0] -Pasily 

Recall  in  the  experiment  that  laminates  of  the  form  [0/90  /0],  n  -  1,2, 3, 4 

n 

and  *  n  "  1,2  wer*  Catted  under  uniaxial  tension  under  room 

temperature  condition.  The  [O^^O^/O^]  laminates  were  also  tested  under 
varying  test  tempsrstures.  In  all  cases,  except  [0/90/0],  transverse  cracks 
in  the  90*-layer  were  recorded  before  final  laminate  failure.  Experimental 
results  obtained  in  this  test  series  have  been  summarized  and  tabulated  in 
Table  3.2,  Section  III. 

According  to  the  finite  element  simulation  procedures  discussed  pre¬ 
viously  for  the  growth  of  transverse  cracking,  energy  release  rate  coefficient 
functions  are  generated  for  each  type  of  laminate  in  this  series.  The 
material  constants  for  the  basic  material  ply  used  in  the  calculation  are 
given  in  Eq.  (3.1)  in  Section  Ill.  Figs.  4.6,  4.7  and  4.8  display  the 
coefficient  functions  C#,  C ^  and  Cy,  respectively  for  the  [0/90n/0]  series. 

In  order  to  use  the  coefficient  functions  for  predicting  the  initiation 
of  transverse  cracking  in  these  laminate,  it  la  necessary  to  define,  at 
the  outset,  the  fracture  quantities  GIc  and  a^. 

As  has  been  already  discussed  earlier,  transverse  cracking  in  the 

90*-layer  is  of  mode-I  fracture.  The  associated  material  resistance  G^c 

must  be  determined  experimentally  in  a  manner  analogous  to  the  cracking 

process  Itself.  The  experiment  discussed  in  Section  3.3(d),  however,  was 

not  appropriate  for  simulating  transverse  cracking.  Rathar,  It  determined 

the  Gyc  value  for  cracking  similar  to  0*/0a  delamination.  In  that  case, 

2 

It  was  found  that  G^O.9  in-lb/in  . 


107 


William  [12],  who  used  Che  AS-3502  graphite-epoxy  system,  found  Chat 

GIe  for  0°/0°  delamination  cracking  was  in  the  range  of  0.8  to  0.92  in-lb/in 

and  for  90*/90*  delamination  cracking  was  in  the  range  of  1.2  to  1.4 
2 

in-lb/in  .  The  latter  simulates  more  closely  with  the  transverse  cracking 
process  discussed  here. 

Since  the  T300/934  material  system  used  in  this  test  series  is  very 

similar  to  the  AS-3502  system  used  by  Williams  [12],  we  shall  use  the  value 
2 

G_  "1.3  in- lb/in  for  all  transverse  cracking  simulations  conducted  in 

1C 

this  study.  Thus, 

Glc  •  1.3  in-lb/in2  (4.21) 

As  for  the  quantity  Sq,  it  is  much  more  difficult  to  determine  experi¬ 
mentally.  Recall  that  the  definition  of  a^  must  be  viewed  at  the  macroscopic 
level  and  in  the. sense  of  the  Griffith's  original  problem  (see  discussion  in 
Section  II).  It  represents  some  effective  flaw  interior  to  the  90°-layer 
with  a  size  of  2aQ. 

Again,  in  the  experiment  discussed  in  Section  3.3(d),  where  specimens 
of  double-side  notches  were  tested,  on  the  basis  of  the  experimental  strength 
of  the  [90g]  specimens,  it  was  suggested  that  the  specimen  had  an  effective 
side  notch  of  size  Sq  ■  2.5  times  the  ply  thickness  (ply  thickness  •  0.0052") 
This  Sq  value  is  associated  with  specimens  having  cut  edges;  man-made  flaws 
along  the  edge  may  have  contributed  greatly  to  the  size  Sq. 
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Figure  4.6  Energy  Coefficient  Function*  C,  for  Transverse 
Cracking  in  [0/90n/0]  Laainataa. 
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Figure  4.7  Energy  Coefficient  Functions  for 

'transverse  Cracking  in  {0/90q/0}  Lanina tea 
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Figure  4.8  Energy  Coefficient  Functions  Cj  for  Transverse 
Cracking  in  [0/90q/0]  Laalnatee 
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In  th«  case  where  the  90*-layer  without  edges  (the  Griffith  problem),  or 
its  edges  ere  constrained  from  cracking,  it  is  believed  that  the  effective 
flaw  should  be  much  smaller.  Thus,  without  a  direct  experimental  determina¬ 
tion,  we  shall  assume  here  that  if  the  thickness  of  the  90* -layer  is  un¬ 
bounded,  an  interior  flaw  of  size  2aQ  “  3t  exists,  which  determines  the  un¬ 
stable  mode-I  crack  propagation  of  the  90*-layer  under  tension.  Or,  alterna¬ 
tively,  the  quantity  a^  associated  with  the  transverse  cracking  simulation  is 
defined  as 

aQ  ■  1.5  t  (4.22) 

t  being  one  ply  thickness  (equal  0.0052"). 

With  GIc  and  a^  defined,  the  prediction  of  the  onset  of  transverse 

cracking  in  this  series  of  tests  are  presented  in  the  following. 

Consider  first  the  laminate  [0/90/0],  n  ■  1.  The  total  thickness  of 

the  90° -layer  is  only  t  -  0.0052".  Any  inherent  flaws  in  the  thickness 

direction  cannot  be  more  than  0.0052"  in  size,  but  can  be  anything  smaller 

than  t.  Thus,  for  all  possible  sizes,  the  largest  possible  energy  release 

rate  is  associated  with  a  *  0.4 1,  see  Fig.  4.6.  The  energy  release  rate 

m 

coefficients  C  ,  C  _  and  C_  at  a  -  a_  are  found  from,  respectively.  Figs, 
e  ei  i  m 

4.6,  4.7,  and  4.8: 

C€(an)  -  136x10%  CeT(affl)  -  39;  CT(aB)  -  2.9xl0-4. 

Then  by  means  of  Eq.  (2.30),  the  maximum  available  energy  release  rate 
for  transverse  cracking  is 

0(«.)  -  IC.UjiJ  +  =.!<*„> S2sT  +  Ct(..)1T21c  <4.23) 
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And,  according  to  the  criterion  for  Initiation,  Eq.  (2.32),  onset  of 


transverse  cracking  Is  defined  by 

6Un)  “  GIc  (4*24) 

Using  the  prescribed  values  of  GJc  *  1.3  in- lb /in2,  t  «  0.0052"  and 
AT  -  TQ  -  Tfcegt  -  300*-75*  -  225*?,  we  calculate  fro®  (4.23)  and  (4.24) 
the  onset  laalnate  strain  e^, 

e%  -  1.03%. 

The  calculated  laalnate  modulus  E  is  14.55  asl.  Hence,  the  onset  laalnate 

x 

stress  o^  Is 

■  150  kai. 

Similarly,  onset  strains  for  the  [0/902/0]  and  [0/ 90^/0]  laalnates 
can  also  be  found  using  the  maximum  strain  energy  release  rate  at  a  ■  afl: 
for  n  •  2,  a  -  0.8t,  and  for  n  •  3,  a  -  1.2t;  all  are  smaller  than 

SB  IB 

Sq  *  1.5t.  Thus,  for  n  *  2 

C#(0.8t)  -  270xl04;  CgI(0.8t)  -  77;  CT(0.8t)  -  5.5xl0“4 
for  n  •  3, 

C#(1.2t)  -  415xl04;  CgT(1.2t)  -  115;  CT(1.2t)  -  8xl0~4. 

The  predicted  onset  strains  are 

*x  "  for  a  "  2 

e^  ■  0.464%  for  n  *  3 

The  computed  laalnate  moduli  Ex  are,  respectively  11.35  and  9.4  xel;  hence, 
the  respective  onset  laminate  stresses  are. 
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a  "  72.6  ksl  for  n  ■  2 
x 


ax  ■  43.6  ksi  for  n  ■  3. 

In  the  ease  of  [0/90^/0] ,  n  ■  4»  the  total  thickness  of  the  90° -layer 
Is  larger  than  2Sq  ■  3t.  Hence,  the  worst  flaw  in  the  thickness  direction 
cannot  be  larger  than  2aQ  -  3t,  Then,  using  the  energy  release  rate 
coefficient  curves  in  Figs.  4.6,  4.7  and  4.8,  we  determine 

Ce(aQ)  -  55Qxl04;  CaT<a0)  -  150;  CT(aQ)  -  O.lxlO-2. 

By  equating  G(aQ)  *  GIc»  we  obtain  the  onset  laminate  strain, 

e  -  0.37Z 
x 

With  the  computed  laminate  modulus  ■  8.12  msl,  the  onset  stress  is 

»  30  ksi. 

The  prediction  for  the  [02/902/02]  laminates  is  carried  out  using  the 
energy  release  rate  coefficients  found  for  the  [0/90/0]  laminate.  But, 
in  this  case,  t  takes  the  value  of  2x0.0052"  for  double  plies  in  each  of 
the  layers.  Hence,  the  predicted  onset  laminate  strain  is 

-  0.63Z 

and  the  corresponding  stress  is 

ax  -  91.6  ksi. 
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As  for  Che  [O^/^^/O^]  laminate,  Che  energy  release  rate  coefficient 
curve  for  [0/902/0]  may  be  used.  Here,  the  length  scale  t  takes  the  value 
2x0.0052"  .  The  value  of  G  Is  evaluated  at  a  •  ag  ■  0.75t  (because  the 
90#-layer  thickness  Is  now  larger  than  aQ) .  Hence,  from  Figs.  4.6,  4.7  and 

4.8,  we  determine  Cg,  and  at  a/t  -  0.75,  from  the  curves  corresponding 
to  n  -  2. 

Ce(aQ)  -  260xl04;  CeT(aQ)  -  70;  CT(a)  -  5.2xl0’4 

and  the  calculated  onset  strain  for  transverse  cracking  is  then 

ex  -  0.3842 

or  the  stress  is 

Sx  ■  43.5  ksi. 

Table  4.1  below  summarizes  the  calculated  results.  Values  in  the 
parentheses  indicate  the  corresponding  experimental  results.  A  comparison 
between  the  calculated  and  the  experimental  results  is  illustrated  in  Fig. 

4.9.  It  is  seen  that  the  transverse  cracking  model  not  only  provided  the 
correct  trend  reflecting  the  effect  of  the  90° -layer  thickness  on  the 
cracking  initiation,  but  also  yielded  reasonably  close  agreement  numerically 
with  the  experiment.  It  should  also  be  noted  that  the  prediction  actually 
gives  the  "minimum"  load  required  for  cracking  initiation,  and  it  is, 
therefore,  the  lower  bound  description. 

In  the  calculations  of  the  onset  laminate  strain  or  stress,  it  is 
noted  that  the  contribution  to  the  energy  release  rate  from  thermal  curing 
stress  is  rather  significant.  Fig.  4.10  shows  for  the  [0/90q/0]  series, 
the  actual  proportions  of  Gg,  GgI  and  G?  in  the  total  energy  release  rate. 
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Ic  Is  seen  that  the  thermal  energy  Itself  is  relatively  small,  especially 
in  the  laminates  of  thin  90°-layer.  On  the  other  hand,  the  thermal/ 
mechanical  coupling  energy  is  the  most  dominant}  it  is  more  pronounced  in 
the  laminates  with  thicker  90° -layer. 


Table  4.1  Summary  of  Calculated  Results 


Laminate 

wmam 

Stress  at  Transverse  Crack 

[0/90/0] 

14.55(13.9) 

150 (No  crack) ,  ksi 

[0/902/0] 

11.35(11.3) 

72.6(59-71) 

[0/903/0] 

9.4(8. 9) 

43. 6  (*+4-53)* 

[0/904/0] 

8. 1(7. 9) 

30  (34-44)* 

[02/902/02] 

14.55(13.9) 

91.6(76-98) 

[02/904/02] 

11.35(11.2) 

43.5(42-57) 

*  Cracks  found  between  the  load  Interval  indicated. 


Finally,  for  the  [02/904/02]  laminates  that  were  tested  at  different 
temperature  environments,  see  Fig.  3.10,  prediction  of  transverse  cracking 
can  be  carried  out  similarly  as  before,  except  now 


AT 


T0" 


T  -  300'F  -  T  . 
test  test 
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(final  failure) 


\ 


[°/90n/0] 


G  *  1.3  in-lb/in 

V 

aQ=  1.5t  =  0.0075* 
At  =  225°F 


Dashed  Lines  are 
Calculated  Results. 


Figure  4.10  Contribution  of  the  Mechanical,  Termal  and  the 
Coupled  Load  Effects  on  the  Total  Energy  Release 
Rate  at  Onset  of  Transverse  Cracking. 


Experiments  [12,  28]  have  suggested  that  the  value  of  Gjc  is  generally 
affected  by  temperature,  that  increases  slightly  with  increasing  temper¬ 
ature  and  vice  versa. 

Lacking  an  experimental  determination  of  the  dependence  of  G  on 

c 

temperature  for  the  T300/934  system,  we  shall  defer  this  question  in  a  latter 

2 

study.  However,  if  we  use  Gjc  ■  1,2  and  G^c  ■  1.4  in-lb/in  ,  respectively 
in  the  calculation,  a  narrow  band  is  obtained  for  the  onset  laminate  stress 
(at  transverse  cracking)  as  a  function  of  the  test  temperature.  Fig.  4.11. 

It  is  seen  that  the  experimental  results  generally  lie  within  the  band, 
except  at  higher  temperature.  The  latter  suggests  that  the  material  GIc 
value  may  actually  be  larger  than  1.4  in-lb/in2  at  temperature  >  120°C. 

4.5  Correlation  with  the  [+25/90n]s-Family 

In  the  experimental  study  on  the  [+25/90  ]  family,  it  is  recalled 

II  3 

that  several  different  cracking  modes  occurred  in  the  laminates.  The  sequence 
of  occurrence  of  the  various  cracking  modes  was  found  to  depend  on  the 
thickness  of  the  90°-layer ,  or  the  parameter  n.  In  order  to  model  these 
crack  growth  events,  we  shall  consider  separately  the  onset  of  transverse 
cracking  and  the  onset  of  delamlnatlon  that  occur  in  the  laminate  series. 

(a)  Prediction  for  Onset  of  Transverse  Cracking.  For  the  [+25/90Q]g  series, 

the  transverse  cracking  mechanisms  are  Identical  to  the  [0/90q/0]  laminates 

that  were  discussed  previously.  The  finite  element  simulation  of  the  crack 

provides  the  energy  release  rate  coefficients  C#,  C#T  and  C^  for  each  value 

of  n(*  1/2, 1,2, 3, 4, 6, 8).  These  coefficient  curves  are  plotted  as  a  function 

of  ths  crack  length  a/t,  t  being  the  thickness  of  a  single  ply  (0.005”), 

Figs.  4.12,  4.13  and  4.14.  It  is  seen  that  these  curves  are  similar  to 

those  shown  in  Figs.  4.6,  4.7  and  4.8  for  the  [0/90/0]  series. 

n 
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Cracking  In  the  [+25/90  ]  Laminates, 
“ns 


Cracking  In  the  1+25/90 ) 


To  predict  the  onset  of  transverse  cracking  for  this  series  of  laminates, 

we  shall  use  the  same  material  constants  as  before.  Namely,  AT  *  225°F, 

2 

Sq  »  1.5t  -  0.0075"  and  G^c  •  1.3  in-lb/in  .  Table  4.2  below  summarizes  the 
appropriate  values  for  C#,  C#T  and  used  in  the  calculation,  and  the 
calculated  onset  laminate  strains  for  transverse  cracking  in  the  series. 


Table  4.2  Summary  of  Calculated  Results  for 
Onset  of  Transverse  Cracking 


[+25/90  1 
ns 

Ce,  *  106 

CeT 

CT,  x  10'4 

V* 

n  -  1/2 

1.0 

41 

2.1 

1.14 

1 

2.5 

75 

6.1 

0.66 

2 

5.4 

140 

10.5 

0.38 

3 

5.7 

161 

10.9 

0.352 

4 

6.0 

154 

9.5 

0.36 

6 

6.4 

148 

8.5 

0.363 

6 

6.5 

140 

7.5 

0.378 

[WSj/’Oj], 

2.45 

73 

6.0 

0.375* 

*  Use  t  -  2  x  0.0052" 


In  the  above  table,  it  is  seen  that  the  thermal /mechanical  coupling 
effect  is  largest  in  the  case  of  n  *  3.  Hence,  the  predicted  onset  laminate 
strain  for  transverse  cracking  is  also  the  lowest.  Also,  in  the  cases  of 
n  -  1/2  and  1,  the  calculated  onset  laminate  atrain  are  large,  because  of 
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chc  available  energy  release  rates  are  small.  In  fact,  for  these  two  cases, 
there  was  no  transverse  crack  found  In  the  experiment  (Instead,  free  edge 
delamination  occurred) .  These  calculated  results  are  graphically  displayed 
in  Fig.  4.15  along  with  the  corresponding  experimental  values  found  in 
Section  III.  It  is  noted  that  the  comparison  between  the  prediction  and  the 
experiment  is  exceptionally  dose,  not  only  in  the  trend,  but  also  in  the 
actual  numerical  value. 

(b)  Prediction  for  Onset  of  Edge  Delamination  It  is  recalled  in  the  experi- 

ment  (Section  III)  that  edge  delamination  was  observed  prior  to  final  failure 

in  each  of  the  laminates  in  the  [+25/90  ]  series.  The  observed  edge  delami- 

—  ns 

nations  generally  occurred  at  loads  in  excess  of  90X  of  the  ultimate 
laminate  strength;  it,  therefore,  may  be  regarded  as  the  initiator  of  the 
laminates'  final  failure. 

In  the  finite  element  modeling,  the  same  ply  properties  were  used 
in  the  calculation,  except  (1)  the  appropriate  critical  energy  release  rate, 
whether  be  it  GIc  or  mixed-mode  G(I,II)c  ;  (2)  the  appropriate  value  for  a^ 
and  (3)  the  ply  properties  of  the  90° -layer  after  it  had  suffered  trans¬ 
verse  cracks. 

As  has  been  discussed  previously,  edge  delamination  along  the  mid-plane 

of  the  laminate  can  be  modeled  as  mode-I  crack.  Experimental  observations 

in  this  type  of  delamlnation,  indicated  that  the  cracked  surface  was  similar 

to  the  crack  surface  of  a  transverse  crack  in  the  90°-layer.  Thus,  for 

2 

mode-I  crack,  the  value  of  G^c  >1.3  in-lb/in  may  be  confidently  used. 


Figure  4.15  Calculated  Laminate  Strain  at  Onset  of  Transferee 
Cracking  in  [+25/90^]^  Laminates. 


On  the  other  hand,  delamination  along  the  25/90  interface  is  generally 


mixed-mode.  As  will  be  shown  later,  the  calculated  to  (Gu ~ 0) 

*  . 

f  ratio  is  approximately  one  to  one  for  laminates  n  »  3  and  n  «  4;  GTTT  is 

F  !  I  lLL 

I  l  much  less  than  G^  for  n  »  2.  Referring  to  the  experimental  data  discussed 
in  Section  111(d)  and  in  particular  the  data  shown  in  Fig.  3.22,  the  possible 
range  of  G^  sbou^  be 


1.5  <  G(I  I1)c  <2.5  in-lb/in 


v4.25) 


The  appropriate  value  for  a^  in  the  free  edge  delamination  cases  should 
be  larger  than  1.5t,  which  was  used  in  the  transverse  cracking  modeling. 

The  reason  is  that  the  crack  is  initiated  from  the  cut  edges  of  the  specimen, 
rather  than  from  the  interior  as  in  the  case  of  transverse  cracking.  In 
the  calculations  that  follow,  aQ  takes  the  value  of  2t, 


a„  -  2t  -  2x0.0052" 


(4.26) 


Finally,  for  the  laminates  of  n  >  2,  multiple  transverse  cracks  were 
formed  in  the  90" -layer  before  edge  delamination.  The  ply  elastic  proper¬ 
ties  of  the  90°~layer  have  obviously  been  degraded.  An  approximate  estimate 
for  the  reduced  90° -layer  stiffness  may  be  determined  by  a  backward  calcula¬ 
tion  using  the  laminated  plate  theory.  In  the  laminate  series  of  n  ,>  2,  the 
following  reduced  property  for  the  90’ -layer  is  used  in  the  calculation: 


-  21.0  Msi 


LT 


LT 


-0.2  ue/*F 


Et  -  1.2  Msi 


TZ 


0.3 


Gtz  -  0.35  Msi 
aT  -  16  ye/°F 


Ez  -  1.72  Msi 


VLZ  *  °’3 


°LZ 


0.94  Msi 


oz  -  16  ve/*F 


(4.27) 
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Note  that  only  four  moduli  are  affected  and  they  are  underlined.  The  global 
effect  of  the  transverse  cracks  gives  approximately  the  same  reduction 
In  the  effective  laminate  stiffness  as  observed  experimentally.  It  will  be 
shown  later  that  the  reduction  of  the  90*-layer  property  is  not  sufficient 
to  predict  the  onset  of  delamination  whenever  transverse  cracking  preceded 
delamlnatlon.  In  these  cases,  the  interactions  between  the  transverse  crack¬ 
ing  and  delamination  processes  become  more  pronounced  as  has  been  discussed 
previously  in  Sections  II  and  III. 

Mid-Plane  Delamination 

Calculation  for  mid-plane  delamination  was  carried  out  for  laminates 
of  n  ■  1/2,  1,2,3  and  4.  The  energy  release  rate  coefficients  associated 
with  the  mechanical,  thermal,  and  the  coupled  components  are  calculated  and 
shown  in  Figs.  4.16,  4.17  and  4.18.  The  general  shape  of  these  coefficient 
functions  .  increase  from  zero  to  a  maximum  value  followed  by  a 

slow  reduction  in  value  as  the  crack  propagates  inward.  This  implies  that 
once  a  mid-plane  delamlnatlon  forms,  it  may  undergo  stable  crack  growth. 

With  the  exception  of  the  case  of  n  “  1/2,  the  magnitude  of  these  shape 
functions  decreases  as  n  is  increased.  In  the  case  of  n  ■  1/2,  the  total 
laminate  is  approaching  the  limiting  case,  say,  of  a  [+25]  laminate,  which 
does  not  delaminate  at  all  in  the  assumed  opening  mode. 

These  energy  release  rate  coefficient  functions  are  calculated  based 
on  the  assumption  that  delamlnatlon  is  the  first  mode  of  damage  in  the 
laminate  and  hence  the  full  ply  properties  are  used  for  all  plies  In  the 
laminate. 
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Figure  4.16  Mechanical  Energy  Release  Rate  Coefficients  for  Mid-Plane 
Delamination  in  the  n  *  1/2,  1,2,3  and  4  Laminates. 


Coupled  Themal/Mechanlcal  Energy  Release  Rate  Coefficients 
for  Mid-Plane  Delamination  in  the  n  ■  1/2,  1,2,3  and  4 
Laminates. 


As  discussed  before,  edge  delamination  occurred  after  the  formation 
of  transverse  cracks  in  the  90° -layer  for  those  laminates  of  n  >  2.  Thus, 
the  post-transverae  cracking  mid-plane  delaminations  are  modeled  a  second 
time  using  the  reduced  effective  moduli  as  given  in  equation  (4.27)  for  the 
cracked  90*-layere  only.  The  corresponding  energy  release  rate  coefficient 
functions  for  n  -  1/2,  1,2,3  and  4  are  shown  in  Figs.  4.19,  4.20  and  4.21  for 
the  mechanical,  thermal,  and  the  coupled  components,  respectively.  These 
coefficient  functions  are  slightly  larger  in  magnitude  than  those  which  do  not 
include  the  effect  of  transverse  cracking.  These  indicate  a  larger  available 
energy  release  rate  for  delamlnatlon.  Thus,  the  inclusion  of  the  effect  of 
transverse  cracks  in  this  manner  will  predict  a  lower  onset  load  for  mid¬ 
plane  delamination. 

For  prediction  of  mid-plane  delamlnatlon,  we  use  AT  -  225°F  and 
GIc  -1.3  in-lb/ln  as  before;  and  Sq  -  2t  as  discussed  earlier.  The  pre¬ 
dicted  onset  laminate  strains  for  n  -  1/2,  1,2,3  and  4  are  summarized  in 
Table  4.3.  It  is  seen  that  the  predicted  onset  strains,  with  or  without  the 
effect  of  reduced  90°-layer  property  (due  to  transverse  cracks)  are  practically 
the  same.  Since,  for  the  cases  of  n  -  1/2  and  1,  the  predicted  onset  strain 
for  delamlnatlon  is  lower  than  the  predicted  onset  strains  for  transverse 
cracking  (see  Table  4.2),  and  since  the  predicted  onset  strains  for  delami¬ 
nation  are  larger  than  the  onset  strains  for  transverse  crack  in  the  cases 
of  n  ,>  2,  the  values  that  are  underlined  in  Table  4.3  govern. 

From  Table  4.3,  it  is  seen  that  the  mechanical  component  of  the  energy 
release  rate  generally  dominates  the  fracture  process  while  the  thermal  energy 
is  quite  insignificant.  Also,  with  the  exception  of  n  -  1/2,  the  onset 
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Calculated  Critical  Strain  and  Energy  Release 
Rates  for  Mid-Plane  Free -Edge  Delaaination. 
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laminate  strain  increases  as  the  thickness  of  the  90*-layer  also  increases. 

This  trend  suggests  that  mid-plane  delamination  may  not  be  possible  when  n  _>  2. 
Indeed,  as  will  be  shown  next,  a  mixed-mode  25/90  de lamination  is  predominant 
when  n  >  2  (this  was  also  observed  in  experiment) . 

Mixed-mode  25/90  Delamination 


The  schematics  of  this  mixed-mode  crack  and  the  finite  element  representa¬ 
tion  have  been  shown  in  Fig.  2.13.  The  energy  release  rate  coefficients  C^, 

CftT  and  CT  for  the  laminates  of  n  •  1/2,  1,2, 3, 4, 6, 8  are  calculated  for  the 
25/90  delamination.  Figs.  22,  23  and  24  show  the  coefficients  calculated 
using  full  ply  properties  (no  consideration  of  transverse  cracking)  for  the 
90#-layer.  And,  Figs.  25,  26  and  27  show  those  calculated  with  reduced  90#-layer 
properties. 

For  predicting  the  onset  strains  for  this  type  of  delamination,  we  used 
all  the  previously  used  material  constants  (AT,  aQ,  etc)  except  that 


G(I,II)c  “  1,5  in“lb/ln  • 


(4.28) 


This  quantity  has  been  suggested  as  the  lower  bound  value  for  mixed-mode 
cracking.  And,  there  is  no  other  experimental  evidence  to  indicate  a  more 
definite  value  for  n)c* 

Accordingly,  the  theoretical  onset  strains  are  calculated;  and  they  are 
tabulated  in  Table  4.4.  Again,  the  effect  of  the  reduced  90°-layer  property 


is  negligible.  Although,  the  predicted  strains  for  n  _>  2  are  uniformly  smaller 
than  that  for  mid-plane  delamination,  the  25/90  delamination  may  not  govern 
when  n  _>  4,  because  the  increasing  trend  of  the  prediction  with  the  Increase 
of  the  90°-layer  thickness  suggests  that  a  different  mode  of  failure  may 
occur  before  this  type  of  delamlnation ,  which  is  due  to  edge  effects  alone. 


bbse?: 


**  ■'f'jjt. fit  ‘-1-&S& 


Figure  A. 23  Coupled  Thermal/Mechanical  Energy  Release  Rate  Coefficients  for 
25/90  Interface  Free-Edge  Delaaination. 


25  Mechanical  Load  Energy  Release  Rate  Coefficients  for  25/90  Interface 
Free-Edge  Delamination  with  Reduced  90°  Ply  Properties. 


Coupled  Thermal/Mechanical  Energy  Release  Rate  Coefficients  Fre e-Edge 
Delaaination  with  Reduced  90s  Ply  Properties. 
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Mixed-Mode  Delamination  From  a  Transverse  Crack  Tip 

From  the  previous  discussions,  it  is  evident  that  formation  of  free  edge 
de lamination  interacts  with  the  formation  of  transverse  cracking  in  the 
1+25/90 i  ]  series.  The  relative  dominance  of  one  over  the  other  depends 

«ea»  U  3 

closely  on  the  parsmeter  n.  For  n  >  2,  transverse  cracks  always  form  before 

delamination.  And,  the  shear  stress  concentration  at  the  transverse  crack 

tip  can  aggravate  delamination  along  the  25/90  interface  where  the  transverse 

crack  is  terminated.  Fig.  4.28  shows  the  interlaminar  shear  stresses  t  as 

yz 

plotted  against  the  distance  from  a  transverse  crack  tip  along  the  25/90  inter 
face.  It  is  seen  that  the  singular  stress  magnitude  near  the  crack  tip  is 
essentially  the  same  for  all  values  of  n.  However,  the  shear  stress  gradient 
away  from  the  crack  tip  becomes  gradual  as  the  value  of  n  increases.  Thus,  in 
terms  of  the  available  strain  energy  for  delamination  along  the  interface, 
the  larger  the  value  of  n,  the  more  available  strain  energy.  Clearly,  for 
a  large  n  (or  thicker  90° -layer),  the  likelihood  of  a  transverse  crack  tip 
induced  delamination  becomes  predominant. 

Fig.  4.29  illustrates  the  schematical  representation  of  the  transverse 
crack  tip  induced  delamination.  Without  considering  other  interacting  factors 
such  a  delamination  growth  can  be  modeled  as  a  1-dimensional  crack  propaga¬ 
tion  in  the  load  direction  along  the  25/90  interface.  Accordingly,  the  finite 
element  simulation,  in  this  case,  is  shown  in  Fig.  4.30.  The  applied  load 
is  represented  by  a  uniform  displacement  6  of  the  +25- layer,  such  that  the 
far-field  axial  laminate  strain  is  a  unity.  Delamination  growth  is  denoted 
by  the  crack  size  a. 


90* -PLIES 


+25* -PLIES 


Figure  4.29  Schematic  of  Transverse  Crack  Tip  Delamination 
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Figure  4.30  Finite  Element  Model  for  Trans¬ 
verse  Crack  Tip  Delamination. 


Numerical  calculations  for  the  energy  release  rate  coefficients  are 
conducted  for  the  laminates  of  n  ■  3,4,6  and  8.  And,  these  are  shown  In 
Fig.  4.31.  These  calculated  coefficients  represent  the  total  sum  of  the 
mode  I  and  11  components;  the  aode-111  component  is  negligible  for  all 
values  of  n. 

It  is  noted  that  for  this  type  of  delamination,  the  modern  shearing 
action  is  predominant.  Initially,  when  the  crack  is  small,  the  crack- tip 
"opening"  is  actually  closed.  As  the  crack  becomes  large  (greater  than  t) , 
then,  the  mating  surfaces  of  the  crack  become  separated  with  some  "opening!*. 
The  emergence  of  the  "opening''  mode  for  each  value  of  n  is  indicated  by  the 
transition  of  the  dashed  lines  to  solid  lines  in  Fig.  4.31. 

Let  us  now  examine  again  the  schematic  representation  in  Fig.  4.29* 

Along  the  free  edge  of  the  [+25/90  ]  laminate,  there  is  still  the  competing 
delamination  action  due  to  the  edge  interlaminar  stresses  on  the  25/90 
interface.  The  available  strain  energy  for  this  edge  crack  action  has  been 
evaluated  earlier  for  the  25/90  mixed-mode  delamination  cases,  and  the 
corresponding  energy  release  rate  coefficients  are  displayed  4n  Figs.  4.25, 
4.26  and  4.27. 

The  driving  forces  (energy  release  rate)  for  the  above  discussed  crack 
axe  mutually  orthogonal  at  the  Instant  of  initiation;  the  transverse  crack  tip 
delamination  is  in  the  load  direction  and  the  edge  del  ami  nation  perpendicular 
to  the  load  direction.  Thus,  in  the  absence  of  a  truly  3-dimensional  crack 
growth  simulation  model,  we  simply  assume  that  the  two  orthogonal  energy  re¬ 
lease  rates  can  be  added  as  a  vector  sum,  which  is  then  the  driving  force 
for  the  triangular  shaped  delamina tion,  as  observed  in  the  experiment  for 


cases  of  n.  >  4 


RELATIVE  CRACK  LENGTH  a/t 


Figure  4.31  Mechanical,  Thermal,  and  Coupled  Energy  Releaee 
Race  Coefficient*  for  Traneverse  Crack  Tip 
Delamination  in  the  n  ■  3,4,6  and  8  Laminates. 


Thus,  the  vector  bub  of  the  energy  release  rate  Is  given  by 


GTOT 


(4.29) 


where  GT(,  is  the  energy  release  rate  associated  with  delanination  emanating 
from  the  transverse  crack  tip,  and  is  the  energy  release  rate  associated 
with  the  mixed -mode  25/90  interface  delaainatlon  at  the  free  edge.  The 
corresponding  initiation  criteria  for  such  a  delamination  is  then  given  by 


Sot^  -G(i,n)c 


(4.30) 


As  in  the  previous  cases  of  mixed  mode  cracking,  we  again  use  values  for 
2 

G^j  xi)c  "  in-lb/in  and  ■  2t.  A  prediction  for  the  onset  of  these 
cracks  are  then  calculated.  .These  calculated  results  are  given  in  Table  4.5 
below. 


Table  4.5  Calculated  Critical  Strain  for  Combined  Transverse 
Crack  Tip  Delamination  and  25/90  Interface  Free- 
Edge  Delamination 


n 

3 

4 

6 

8 

«> 

0.611 

0.596 

0.573 

0.555 

Figure  4.32  summarises  the  predicted  onset  strains  of  edge  delamination 
for  the  series  of  laminates,  along  with  the  experimentally  obtained  values. 
The  different  symbols  used  in  the  figure  distinguish  the  various  delsmlnation 
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Figure  4.32  Comparison  of  Calculated  Results  with  Experiment,  [+25/90^1  Series 


modes.  It  aiou  that  the  analysis  captures  the  physical  mechanisms  quite 
accuratelyt  although  the  quantitative  prediction  showed  some  discrepancy 
with  the  experiment,  especially  for  n  >  4,  In  the  latter  case*  It  is  believed 
that  the  reason  for  the  low  laminate  strain  at  onset  of  delamination  is  due 
to  the  excess  strain  energy  released  during  the  unstable  propagation  of  the 
transverse  cracks  (see  Fig.  2.8  and  the  associated  discussion  in  Section  II). 
This  excess  energy  in  laminates  of  n  ■  6  and  8  is  very  large;  the  associated 
dynamic  effect  is  believed  to  cause  the  prematured  delamlnatlon  growth,  which  la 
emanating  from  the  transverse  crack  tip.  Indeed,  if  we  examine  the  experi¬ 
mental  data  shown  in  Fig.  3.11,  we  see  for  example  that  in  the  laminate  of 
n  *  8,  delamlnatlon  along  the  25/90  interface  followed  Immediately  the  for¬ 
mation  of  transverse  cracking,  while  the  final  laminate  failure  followed 
immediately  the  delamlnatlon.  It  seems  that  the  sequence  of  events  from 
transverse  cracking  to  final  failure  is  almost  continuous  whenever,  the 
90°-layer  is  excessively  thick. 

4.6  Mixed-Mode  Fracture  in  [Os] -Family  . 

In  this  sub-section,  we  shall  present  some  of  the  numerical  details  in 
simulating  the  crack  growth  emanating  from  the  notch-tip,  along  the  fiber 
direction  for  the  [0g]  family.  In  this  case,  the  same  finite  element  pro¬ 
cedures  were  followed,  except  that  the  basic  formulation  of  the  finite 
element  stiffness  matrix  was  based  on  a  "plane  stress"  approach,  rather  than 
the  "generalized  plane  strain"  approach.  Thus,  the  mode  of  cracking  generally 
consists  of  mode  I  and  moda  II  only.  The  ratio  between  and  depends 
on  the  fiber  angle  8. 
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Fig.  4.33  shows  Chs  calculated  energy  release  rates  Gg  and  G^  as  plotted 

against  the  crack  growth  size  a  (see  Fig.  3.21  for  the  crack  growth  geometry) . 

2 

The  unit  of  G  shown  in  the  figure  is  in-lb/in  per  unit  far-f laid  applied 
strain. 

From  Fig.  4.33,  it  is  seen  that  node  1  action  is  predominant  for  fiber 
angle  9  >  45°;  and  the  mode-II  component  of  G  is  practically  negligible. 
G-curves  for  9  >  45*  are  monotonlcally  increasing  with  the  crack  size  a, 
suggesting  an  unstable  crack  growth  behavior.  As  the  fiber  angle  decreases 
from  45*,  mode-II  contribution  increases  rapidly.  In  particular,  when  9  <  30* 
the  energy  release  rate  curves  show  a  decrease  as  the  crack  grows  from  the 
notch  tip.  As  the  crack  grows  further,  the  available  energy  then  increases 
slightly.  Thus,  a  minimum  value  of  G  exists  for  0.01"  <  a  <  0.04".  This 
behavior  of  the  energy  release  rate  Indicates  a  stable  crack  growth  for  small 
fiber  angles,  say  9  <  20°. 

Assuming  that  a  macroscopic  flaw  of  size  about  0.01"  exists  near  the 
notch  tip  along  the  fiber  direction  (this  is  the  same  as  a0«v2t  as  before) , 
then,  the  values  of  Gg  and  taken  at  a  *  &q  will  determine  the  onset  of 
crack  growth  along  the  9-direction.  These  values  of  Gj.  and  GII  were  shown 
earlier  in  Fig.  3.20;  and  they  were  used  to  correlate  with  the  experimental 
data  (Fig.  3.19).  In  the  process,  we  determined  the  material  resistance 
G(i  under  the  varying  degree  of  mixed-mode  crack  actions.  As  la  shown 
in  Fig.  3.22,  gg^c  shows  a  dependence  on  the  ratio  when  it  is 

measured  at  the  macroscopic  level. 
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Components  of  Double  Side-Notched  -  [9gJ-Family 


V.  CONCLUDING  REMARKS 


In  this  report,  we  have  presented  the  research  results  obtained  in  the 
AFOSR  sponsored  program  under  the  contract  F49620-79-C-0206 . 

Before  this  report,  a  total  of  six  (6)  technical  papers  and  a  Ph.D. 
dissertation  were  published  during  the  contract  period  from  1  September  1979 
to  30,  September  1981.  These  papers  have  dealt  individually  some  of  the 
major  developments  accomplished  in  this  program,  see  Ref.  [29-35]. 

At  the  conclusion  of  this  report,  the  following  observations  are  made. 

(a)  Sub-laminate  crack  growth,  such  as  transverse  cracks  and  edge  delamina¬ 
tion  discussed  in  this  report,  is  basically  a  fracture  failure  process.  It 
has  been  demonstrated  that  the  phenomenological  description,  via  ply  elasticity 
and  the  classical  fracture  mechanics,  can  be  applied  successfully  to  model 
some  of  the  often  observed  cracking  processes. 

But,  it  must  be  pointed  out  that,  given  a  laminate  under  load,  there 
may  be  a  multitude  of  crack  modes  that  interact  with  one  another.  The  result¬ 
ing  complex  network  of  crack  propagation  may  present  a  formidable  mathematical 
problem  if  it  is  also  to  be  described  by  the  same  fracture  analysis  approach. 

Rather,  we  have  learned  from  the  simpler  analysis  the  important  material 
and  geometrical  parameters  which  have  profound  influences  on  the  sub-laminate 
crack  growth  behavior.  For  example,  the  effect  of  the  90°-layer  thickness  on 
the  transverse  cracking  and  delamlnatlon  processes  deserves  due  attention 
in  practical  design.  Similarly,  the  material  toughness  Gc  is  another  param¬ 
eter  which  should  be  considered  in  composite  laminate  design. 
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(b)  Sub-laminate  crack  growth  generally  lowers  the  strength  of  the  laminate. 
The  effect  pertains  to  the  growth  stability  of  the  sub-laminate  cracks.  In 
the  [±25/90n]a  series  studied  in  this  report,  see  Fig.  3.11,  it  la  seen  that 
the  failure  strain  of  the  load-bearing  ply  (the  +25° -layer)  decreases 
dramatically  when  n  >_ 4.  The  in-sltu  stresses  in  the  +25° -plies  just  before 
failure  varies  profoundly  with  the  value  of  n,  as  shown  in  Fig.  5,1.  Clearly, 
the  failure  process  in  the  load-carrying  +25°-plies  is  again  a  fracture  event. 
This,  however,  has  not  been  analyzed  in  this  report. 

(c)  The  flnlte-element/crack-closure  procedure  is  a  useful  computational 
simulation  technique.  In  view  of  the  dimensional  requirement  of  the  phenomeno¬ 
logical  model,  it  is  felt  that  a  simple  formulation  (without  inclusion  of 
stress  singularity,  for  instance)  is  sufficient  for  purpose  of  simulating  the 
kinds  of  crack  propagation  such  as  discussed  in  this  report. 

(d)  It  is  felt  that  the  determination  of  G  associated  with  sub-laminate 

c 

cracks  has  been  Inadequate.  In  particular,  the  process  of  mixed-mode  failure 
remains  almost  unexplored.  Continued  and  future  researches  may  provide  more 
information  soon. 
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